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Abstract 

All possible orbital trajectories and their analytical expressions in the 
Schwarzschild metric are presented in a single complete map characterized 
by two dimensionless parameters. While three possible pairs of parameters 
with different advantages are described, the parameter space that gives the 
most convenient reduction to the Newtonian case is singled out and used 
which leads to a new insight on Newtonian limits among other results. 
Numerous analytic relations are presented. A comparison is made with 
the widely used formulation and presentation given by S. Chandrasekhar. 

PACS numbers: 04.20.Jb, 02.90.+p 



1 Introduction 



As Einstein's theory of general relativity is approaching its centennial celebra- 
tion, the problem of understanding the orbital trajectories of a particle around a 
very massive object still retains great interest. Although current research has fo- 
cused on the properties of spinning black holes and on the effects of gravitational 
waves, the original analytic solutions for the Schwarzschild metric are still used 
for the study of astronomical objects for which the assumption of static spher- 
ical symmetry is reasonable. The orbital solutions to the Schwarzschild metric 
also provide useful guides from which various perturbation methods originate. 

Although the orbital problem for the Schwarzschild metric has been solved 
in terms of the equations for the orbits, a complete picture of the solution space 
in terms of a single clear map has been lacking. Given, say, the trajectory 
of a particle around some star or black hole, one would be interested in the 
characterizing parameters of such a trajectory and how it relates to other pos- 
sible trajectories. Indeed, the question of constructing a universal map and 
what physical parameters should be used for the coordinates (the "latitude" 



and " longitude" ) of such a map did not seem to have been seriously considered 
before our recent work [1-4]. It is well known [see e.g. ref.5], for example, that 
the discriminant of a relevant cubic equation is used to establish the criteria 
for distinguishing different types of possible orbits in the Schwarzschild metric. 
However, because of the absence of a universal map, clear boundary curves that 
separate different regions of the parameter space for different types of orbits 
and trajectories were never clearly presented before our work in ref.4. 

The most widely used analytic solutions and analyses for orbits and tra- 
jectories in the Schwarzschild metric were based on Chandrasekhar's work [5] 
that originated with Darwin [6]. We shall make a comparison of our analysis 
with theirs at the end of this paper. Two notable features of Chandrasekhar's 
analysis are that (i) his analytic solutions for the orbits are implicit, not ex- 
plicit, expressions that include a parameter that is related to the polar angle 
in a somewhat complicated implicit way, and (ii) one of his two parameters for 
characterizing the orbital trajectories can be real or purely imaginary, making 
the regions for different types of trajectories difficult to be placed on a map. 
Among the earlier work, we should mention the publications of Forsyth [7] and 
Whittaker [8] who gave some special cases of the analytic solutions for the orbits 
that we present here (and in refs. [1-4]) but in different forms. Hagihara [9] 
had a detailed analysis of all possible orbits in terms of Weierstrassian elliptic 
functions but his analysis involved some lengthy algebraic expressions and his 
division of all possible orbits into seventeen cases also unnecessarily complicates 
the picture. But he did present a map with two dimcnsionlcss parameters as 
coordinates [9, Fig. 13] even though he failed to clearly identify the physical sig- 
nificance of the two parameters, and he classified the orbits largely according to 
the roots of the relevant cubic equations. 

Two important conclusions to be drawn from our derivation of analytic solu- 
tions for the orbits and trajectories in terms of the Jacobian elliptic functions are 
the following: (a) The consequence of having an analytic solution for the orbits 
is not just that it is exact and convenient, as a numerical solution of the differ- 
ential equation could easily produce the same precise trajectory. The more im- 
portant consequence is that it cleverly picks out the appropriate characterizing 
dimensionless parameters for the orbits and these parameters are dimensionless 
in terms of some combinations of physical quantities that could not have been 
easily guessed otherwise; and (b) The well-studied Jacobian elliptic functions 
allow many simple relations among different parameters for a number of special 
cases that are very useful for checking numerical computations. It is with these 
considerations in mind that we present the characterization of all possible orbits 
and trajectories on a map with the special choice of two dimensionless physical 
parameters as the coordinate axes. 

In this paper, we summarize and extend the most essential features and 
expressions of our previous analysis and present them in a complete picture. 
We describe three pairs of dimcnsionlcss parameters that can be used for the 
coordinates of the map of all possible orbits and trajectories, each of which 
exhibits some special advantages. Among the useful outcomes are the following: 
(1) All possible orbital trajectories of a particle in the Schwarzschild metric are 
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placed on a universal map; (2) The boundary between the two regions that 
separate different types of trajectories is clearly shown on the map; (3) One 
particular set of parameters serves to illuminate some interesting properties 
with regard to the meaning and existence of Newtonian limiting orbits. We also 
present an abundance of simple analytic results and relations for some special 
cases in various parts of the map that are useful for many purposes, not the 
least of which is to provide a way to check results from numerical computations. 

Here is a brief outline of this paper. In Section 2, we introduce the fundamen- 
tal differential equations in the Schwarzschild geometry and three possible pairs 
of characterizing parameters. In Section 3, we present three explicit analytic 
solutions for all possible orbital trajectories of a particle in the Schwarzschild 
metric, define two regions and give examples of orbital trajectories in the two 
regions. While these three fundamental solutions for the orbits may have been 
given in different forms by other authors, we have explicit analytic expressions 
for many related relevant physical quantities, especially those for orbits of the 
hyperbolic type, that were not given by other authors. One of these expressions 
(for the impact parameter of a particle) corrects an erroneous expression given 
by Chandrasekhar [5]. In Section 4 we describe the boundaries of the two re- 
gions in the forms of curves expressed by simple equations. In Section 5, we 
present analytic expressions for the orbital trajectories on these boundaries and 
provide examples. In Section 6, we present a clear picture of how the orbital 
trajectories change across the boundaries of the two regions. In Section 7, we 
introduce a large collection of analytic relations that relate various physical pa- 
rameters at various special locations on our map. These have been very useful 
for checking our numerical calculations and can be useful for checking other nu- 
merical computations. However, this section could be skipped on a first reading. 
In Section 8, we present the Newtonian and non-Newtonian correspondences of 
our general relativistic results and show a sector of a region on our map in 
which there exist elliptical orbits that have properties which have not seemed 
to have been previously recognized. We also derive an approximate formula for 
the general relativistic correction to a Newtonian hyperbolic orbit. Our calcu- 
lation gives a slightly different result from that given by other authors for a 
proposed experiment involving a spaceship that can be carried out in our solar 
system. In Section 9, we relate the parameters that we use in our analysis with 
those used by Chandrasekhar [5] . This comparison can be used for reference by 
those readers who are familiar with Chandrasekhar's work. In Section 10, we 
summarize the principal results of this paper. 

2 The Differential Equations and the Charac- 
terizing Parameters 

We consider the Schwarzschild geometry, i.e. the static spherically symmetric 
gravitational field in the empty space surrounding some massive spherical object 
such as a star or a black hole of mass M. The Schwarzschild metric for the empty 
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spacetime outside a spherical body in the spherical coordinates r, 0, (f> is [10] 



dl 2 =c 2 (l-^dt 2 - (l-") dr 2 -rW-r 2 sin 2 M/> 2 (1) 



where 



2GM 

a =—^~ ( 2 ) 



is the Schwarzschild radius, G is the universal gravitation constant, and c 
is the speed of light. If [x M ] = (t,r,9,(f>), then the worldline x' x (t), where r is 
the proper time along the path, of a particle moving in the equatorial plane 
= 7r/2, satisfies the equations [10] 

(i-")i = «, (3) 



r 2 <f> = h, (5) 

where the derivative ■ represents d/dr. The coordinates r and (/> describe the 
position of the particle relative to the star or black hole situated at the origin. 
The constant h is identified as the angular momentum per unit rest mass of the 
particle, and the constant k is identified to be the total energy per unit rest 
energy of the particle 

k=— 2 , (6) 
mocr 

where E is the total energy of the particle in its orbit and m is the rest 
mass of the particle at r = oo. Substituting eqs.(3) and (5) into (4) gives the 
'combined' energy equation 

,2 . „2 



.2 h 
r +C + ^2 v 



( 1 - a ) - ^ = c 2 « 2 . (7) 



Substituting dr/dr = {dr / dcj)){d<j) / dr) = (h/r 2 )(dr /d<j>) into the combined 
energy equation gives the differential equation for the orbit of the planet 



du^ 2 
dtp 



an* ~u 2 + (2GM/h 2 )u + c 2 { K 2 -l)/h 2 (8) 



where u = 1/r. We define the dimensionless distance q of the particle from 
the star or black hole, measured in units of the Schwarzschild radius, by 



r 1 I 

q= - = — = -. 9 
a au v 



In terms of v, eq.(8) becomes 



4 



dv \ 2 



d<j)J 



= v 6 -v 2 + 4s 2 v + 4s 2 {n 2 - 1), (10) 



where 
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By changing the variable from v to another dimcnsionless quantity U dchned 

by 

1 (a 1\ If 1\ 1/1 1 



A\r 3 J 4V 3/ A\q 3 
or 

v=\ = \+W, (12) 

eq.(10) becomes 

/rlTJ\ 2 

= AU i -g 2 U-g 3 (13) 



dU'- 2 



where 



1 2 

■ 92= 12- S 

and where we have defined the dimensionless quantity e 2 given by 

6 - 1+ (GM) 2 - 1+ S 2 ■ I 15 ) 

It is clear that as soon as we use the dimensionless distance q (and the 
dimensionless angle (f>) to express the differential equation (8) in the form of 
eq.(10) or (13), the solutions for the orbits should be conveniently characterized 
by a pair of dimensionless parameters (k 2 ,s 2 ) or (g2,gs)- Indeed, g 2 and g 3 
are called the invariants of the Weierstrassian elliptic functions [11] that would 
represent the solution for U for the differential equation (13). Notice that while 
the constant of motion k defined by eq.(6), which is the total energy per unit 
rest energy, can obviously be one of the dimensionless parameters, the other 
constant of motion h which is the angular momentum per unit rest mass is not 
a dimensionless quantity and should not be used by itself as a parameter for 
representing the orbits in a universal map. That s defined by eq.(ll) could 
be used as the second characterizing parameter was recognized by one of us in 
ref.[l]. In fact, the orbit solutions depend on k 2 and s 2 as shown by eqs.(10), 
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(13) and (14). While we shall briefly discuss the parameter spaces described 
by the coordinates (k 2 , s 2 ) and (32,53), for most of this paper we shall use the 
parameter space characterized by the coordinates (e 2 , s 2 ). The principal reason 
for this choice of parameters is that the parameter e conveniently reduces to the 
well known parameter called eccentricity for the classical Newtonian case. In our 
earlier work [1,2], there was an oversight on our part in which we unnecessarily 
limited ourselves to the case e 2 > 0, but we remedied it in ref.4. 

We call the two dimensionless parameters e 2 and s 2 defined by eqs.(15) and 
(11) the energy eccentricity and gravitational field parameters, or simply the 
energy and field parameters, respectively. They are the principal parameters 
we use as coordinates for characterizing the orbital trajectories of a particle in 
the gravitational field produced by the massive object or black hole of mass M. 
We use the parameter space (e 2 ,s 2 ), and (e 2 ,s) as well, where —00 < e 2 < 00, 
s > for classifying and identifying all possible types of orbital trajectories. 
In doing so, we shall put all possible trajectories in the Schwarzschild metric 
on a two-dimensional universal map on which the two coordinate axes are real 
and represent dimensionless physical quantities. For e 2 > 0, the three cases of 
0<e<l,e = l and e > 1 can be used to classify the orbits as elliptic- type, 
parabolic-type and hyperbolic-type, similar to those in the Newtonian cases; For 
e 2 < on the other hand, the orbits do not have any Newtonian correspondences 
and this will be explained and discussed in Section 8. 

3 Three Analytic Solutions for the Orbital Tra- 
jectories in Two Regions 

We consider eq.(13) and derive its analytic solutions in this section. 
The discriminant A of the cubic equation 

W 3 -g 2 U~g 3 = (16) 

is defined by 

A = 27g 2 - g 3 2 . (17) 

For the case A < 0, the three roots of the cubic equation (16) are all real. 
We call the three roots ei,e 2 ,e 3 and arrange them so that e\ > e 2 > e 3 ; the 
special cases when two or three of the roots are equal will be considered also. 
For the case A > 0, the cubic equation (16) has one real root and two roots that 
are complex conjugates. The analytic solutions of eq.(13) that we present below 
will give the distance r of the particle from the star or black hole in terms of the 
Jacobian elliptic functions [11] that have the polar angle <j) in their argument 
and that are associated with a modulus k that will be defined. 

There are three relevant analytic solutions of eq.(13). The first two solutions 
are for the case A < 0, and the third solution is for the case A > 0. The region 
that contains all the coordinate points (e 2 ,s 2 ) for the first two solutions for 
the case A < is called Region I (seefTl), and the region that contains all the 
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coordinate points (e 2 ,s 2 ) for the third solution for the case A > is called 
Region II. There are two physically inadmissible regions: the region below the 
horizontal axis (s 2 < 0) and the region corresponding to k 2 < which is the 
region (called Region II' in Fig.l) above the curve s 2 > s\ where s 2 will be 
given in the next section. The boundary between Regions I and II consists of 
two segments that we call s 2 = s\ (the upper boundary of Region I) and s 2 = s' 2 
(the left boundary of Region I) that will be defined in the next section. The 
three boundary curves s 2 = 0, s 2 = sf and s 2 — s\ all satisfy A = 0. 

We now present the three relevant analytic solutions of eq.(13). 

Solution (Al) For A < 0, e\ > e 2 > U > e 3 (Applicable in Region I). 

Writing the right-hand side of eq.(13) as 4(ei — U)(e 2 — U)(U — e 3 ), the 
equation for the orbit is 



- = - + 4e 3 +4(e 2 - e 3 )sn 2 ( 7 0, k) 
q o 

1 . / .1- cn(2j(b, k) 

= 3 + ^ + ^- e * \ + M2^ky ^ 

< 1, the initial point at cf) — gives U — e 3 (i.e. 1/q = 
1/3 + 4e 3 ) and thus dU/dcj) = from eq.(13) which means that the particle at 
this point is in the direction perpendicular to the radius vector. For the case 
e 2 > 1, the initial point is discussed in eq.(28). The constant 7 appearing in 
the argument, and the modulus k, of the Jacobian elliptic functions are given 
in terms of the three roots of the cubic equation (16) by 



7=(ei-e 3 ) 1/2 , (19) 
&i - e 3 
d - e 3 ' 



fc 2 = ^^. (20) 



where ei, e 3 are given by 



and where 



1/2 

d = 'l I — I cos 



e3 



(S) 



3 

, ;.\l/2 (0 2- \ 



2?x 1/2 



cos = g 3 ^-3 J . (22) 

The modulus k of the elliptic functions has a range < k 2 < 1. For k 2 = 0, 
sn(x, 0) = sinx, cn(x, 0) = cosx, dn(x, 0) = 1. For fc 2 = 1, sn(70, 1) = 
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tanh(70), cn(j(f), 1) = dn(j(j),l) — sec h(-f(f)). The period of cn{2 r )(j), k) isAK(k), 
and the period of dn(2-f(f>,k) and of sn 2 (7^, k) is 2K(k), where K(k) is the 
complete elliptic integral of the first kind [11]. As k 2 increases from to 1, 
K(k) increases from tt/2 to 00, and k 2 is a useful parameter for the orbit. 

The part showing Region I in [T] is enlarged in [2] in which the left boundary 
of Region I is denoted in [l] by the curve s[ for which k 2 = 0, and the upper 
boundary of Region I is denoted in[l]by the curve Si for which k 2 = 1. The 
intersection point V of s[ and si, which we shall call the vertex, is (e 2 ,s) = 
(—1/3, l/\/l2) and corresponds to ex — e 2 = e 3 = 0. It is a special point of 
interest. Three other constant k 2 curves for k 2 = 0.01, 0.2 and 0.5 in Region I 
are also shown in [2] The constant k 2 curves in Region II will be discussed later. 

A typical orbit given by eq.(18) (not on any one of the three boundaries) in 
Region I is a precessional elliptic- type orbit for e 2 < 1, a parabolic- type orbit for 
e 2 = 1, and a hyperbolic-type orbit for e 2 > 1 [2,3]. For the elliptic-type orbits 
(e 2 < 1), e 3 > —1/12, the maximum distance r max (the aphelion) of the particle 
from the star or black hole and the minimum distance r m j n (the perihelion) of 
the particle from the star or black hole, or their corresponding dimensionless 
forms (? max (= f max /a) and q m i n (= r m ; n /a), are obtained from eq.(18) when 
70 = and when 7</> = K(k) respectively, and they are given by 

^+4e 3 , (23) 



Qmax 3 

and 

Z*- = 1 + ^ (24) 

ymin *-> 

where e 2 and e 3 are determined from eqs.(21), (22) and (14) in terms of e 
and s. The geometric eccentricity e of the orbit is defined in the range < e < 1 
by 

(25) 



^max "T ^*min Qmax "T ^min 

If we use (j max and q min given by eqs.(23) and (24) and write 

£= l/6 + e 2 + e 3 ' (26) 

then the definition of e can be extended to e > 1 but it should be emphasized 
that any e outside the range between and 1 does not have the physical meaning 
given by eq.(25). We shall have more to say about this point in Section 9. 

The precessional angle A</> is given by 

W= 2 JW-2*. (27) 

7 

Whether the orbit is closed (and hence periodic) or non-closed (and hence 
non-periodic) depends on whether A0/7T is or is not a rational number. 
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An example of a precessional elliptical-type orbit for e 2 = 0.25 and s = 0.2 
in Region I is shown in [3] 

For the parabolic-type orbits (e 2 = 1), e 3 = —1/12, and g max = oo. The 
orbits are unbounded. The particle is initially at infinity at <f> — and initially 
moves counterclockwise in a direction perpendicular to the line joining it to the 
origin. The incoming trajectory coming from infinity at <f> = goes around 
the origin counter-clockwise according to eq.(18) as <p increases from and the 
outgoing trajectory makes an angle <j) — 2K(k)/j with the horizontal axis in 
going to infinity. A particle coming from infinity can go around the star or black 
hole many times as its polar angle <f> increases from to 2K{k)/ r ) before going 
off to infinity. 

An example of a precessional parabolic-type orbit for e 2 = 1 and s = 0.4 in 
Region I is shown in[4j 

For the hyperbolic-type orbits (e 2 > 1) [2,3], e-$ < —1/12, the initial point 
of the particle is not at cf> — but is at infinity at a polar angle <f> = ^i, where 

is given by 



sn(7\E'i, k) 

or 

and where 7 and k are defined by eqs.(19) and (20) 

2K(k) 
7 

is the next value of <f> for q to become infinite. Thus eq.(18) gives the 
hyperbolic-type orbit in Region I for 

*1<0<*2- (31) 

The orbit given by eq.(18) for e 2 > 1 describes a particle approaching the 
star from infinity along an incoming asymptote at an angle ^1 to the horizontal 
axis, turning counter-clockwise about the star to its right on the horizontal axis, 
and leaving along an outgoing asymptote at an angle ^2 ■ A particle approaching 
along the incoming asymptote can go around the star or black hole many times 
as its polar angle 4> increases from ^1 to ^2 before going off to infinity along the 
outgoing asymptote. The x coordinates of the points where the two asymptotes 
intersect the horizontal axis relative to the origin are given by 




(28) 



(29) 

As 6 increases from 
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1 1 

Xal ~ ~2s 2 v / ^Tsin*i' 
11 

Xa2 ~ 2s*V^l^*~2' (32) 

A negative (positive) sign for the x a i means that the corresponding asymp- 
tote intersects the horizontal axis to the left (right) of the origin. The minimum 
distance g m i n of the particle from the star or black hole is obtained from eq.(24) 
when the particle is at a polar angle 4> = K(k)/j, where dr/d<f) = at q = q nlln . 
A straight line through the origin O (where the star or black hole is located) 
that makes an angle 

» - ™ (33) 

7 

with the horizontal axis is the symmetry axis of the hyperbolic-type orbit. 
Let b be the impact parameter, i.e. the perpendicular distance from the origin 
(where the star or black hole is located) to the asymptote. The dimensionless 
impact parameter b/a, where a is the Schwarzschild radius, can be expressed 
as 

- = r (34) 

An example of a precessional hyperbolic- type orbit for e 2 = 4 and s = 0.2 
in Region I is shown in [5] 

The orbital trajectories that occur on the left and upper boundaries of Region 
I are special cases that deserve separate attention and will be discussed in Section 
5. 

We now present the second solution. 

Solution (A2) For A < 0, U > e\ > e% > es (Applicable in Region I). 
We write the right-hand side of eq.(13) as 4(U - ei)(U - e 2 )(U - e 3 ). The 
equation for the orbital trajectory is 

1 1 ei - e 2 sn 2 (j<j), k) 

" = q+ 4 2i a u\ ' ( 35 

q 6 an ,yy<p,k) 

where 7, k, ei, e 2 and are given by eqs.(19)-(21) as in the Solution (Al). 
This solution gives a terminating orbit. The point at 4> = has been chosen to 
be given by 

- = \+le x . (36) 
qi 6 

The particle, starting from the polar angle <f> — at a distance q\ from the 
black hole, plunges into the center of the black hole [2] when its polar angle <j>i 
is given by cn(7</>i, k) — 0, i.e. when 
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K(k) 



(37) 



i = 



7 



Considering Solutions (Al) and (A2) together, each point (e 2 ,s) of the 
parameter space in Region I allows two distinct orbits, one of which is non- 
terminating (that consists of three types: elliptic, parabolic and hyperbolic) 
and the other terminating. At the same coordinate point, quantities that de- 
scribe the two distinct orbits are related. For example, by noting e± + e2 + e% — 
and from eqs.(23), (24) and (36), q\ can be expressed as 



where q m i n and q max are the minimum and maximum distances for the non- 
terminating orbit at the same coordinate point (e 2 , s 2 ). It will be noted that q\ 
is less than q m im i-e. for the terminating orbit the particle is assumed initially 
to be closer to the black hole than the <7 m i n for the associated non-terminating 
orbit, except at k 2 = 1 where q\ = q m ; n and the particle has a circular instead 
of a terminating orbit that will be explained later. We note that q\ is finite even 
for e 2 > 1. 

Three examples of these terminating orbits in Region I for the same coordi- 
nate points of s = 0.2 and e 2 = 0.25, 1, and 4 are shown in |6j [7] and [8] (These 
were also the coordinate points used for [3j [4] and [5| . Notice that the initial dis- 
tances of the particle from the black hole given by qi are all finite and are less 
than <7 m i n for the non-terminating orbits. This corrects an erroneous statement 
made by us in ref. 3 after eq.(75). 

We now present the third solution. 

Solution (B) For A > (Applicable in Region II). 

Define 



where 53 and A are defined by eqs.(14) and (17). The real root of the cubic 
equation (16) is given by 




(38) 




(39) 



a = A + B (40) 

and the two complex conjugate roots b and b are — (A+B)/2±(A — B)\fiij2. 
We further define 



7= [3(A 2 + AB + B 2 )) 1 ' i 



(41) 
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and 



fc2 = 1 _ 3(A + B) = 1 _ 3a 

2 4^3(A 2 +AB + B 2 ) 2 4 7 2 ' 

The curves of constant fc 2 , the squared modulus of the elliptic functions 
that describe the orbital trajectories, radiate from the vertex point V whose 
coordinates are (e 2 , s 2 ) = (—1/3, 1/12) (see[2|. In Region I, fc 2 increases counter- 
clockwise as one goes from the curve k 2 = for s — s[ to the curve k 2 = 1 for 
s = s l5 and as one keeps going counter-clockwise to Region II, k 2 decreases from 
k 2 = 1 and returns in value to k 2 — 0. The squared modulus k 2 is obtained 
from eq.(20) for Region I and from eq.(42) for Region II. Two k 2 = const, curves 
in Region II are rather special: k 2 = 1/2 + \/3/4 = 0.933012 is a horizontal line 
to the right of V and k 2 = 1/2-y/S/A = 0.0669872 is a horizontal line to the left 
of V. The k 2 = 1/2 curve in Regions I and II has a simple analytic expression 
that will be given later. 

Writing the right-hand side of eq.(13), with U > a, as A(U-a)(U-b)(U-b), 
the equation for the orbital trajectory is 



1 1 , A 2 l-crc,(27<£,fc) 

- = - + 4a + Arc — 

q 3 7 1 + cn(2j4>, k) 

= 1 + 4a + 47 2 in 2 (70, k)dn 2 (jcf), fc). (43) 
3 

This is a terminating orbit and thus Region II allows only terminating orbits. 
For e 2 < 1, a is greater than —1/12, the initial distance of the particle at 
4> = (in a direction dq/d<p — 0) is given by 

1 = I+4a. (44) 

02 3 

The particle plunges into the center of the black hole [2,3] when its polar 
angle (f> = (f> 2 is given by 

* - (45) 

7 

where 7 and fc are given by eqs.(41) and (42). 

An example of a terminating orbit for e 2 = 0.25 and s = 0.3 in Region II is 
shown in [9] 

For e 2 = 1, a = —1/12, and the initial distance of the particle at cj> = is 
at infinity (in a direction dq/d<p = 0). An example of a terminating orbit for 



e 2 = 1 and s = 0.3 in Region II is shown in 10 in which the initial position of 
the particle coming from infinity is at = 0. 

For e 2 > 1, a is less than —1/12, the particle comes from infinity at the polar 
angle <f> = $2, where $2 is given by 

cn(2 7 $ 2 ,fc) = (46) 
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where 



and where 7 and k are given by eqs.(41) and (42). The orbit given by 
eq.(43) describes a particle coming from infinity along an incoming asymptote 
at an angle $2 to the horizontal axis, and as 4> increases, the orbit terminates 
at the center of the black hole [3] at a polar angle 2 given by eq.(45). There is 
thus only one asymptote, and the x coordinate of the point where the asymptote 
intersects the horizontal axis relative to the origin is given by 

(47) 



2s 2 v 7 ^^! sin $ 2 ' 

and the dimensionless impact parameter is given by eq.(34). 

An example of a terminating orbit for e 2 = 4 and s — 0.3 in Region II 
is shown in [TT] in which the particle comes from infinity along the asymptote 
shown by the dotted line. 

We have given three examples of these terminating orbits in Region II which 
are shown in [9] [To] and for the coordinate points for e 2 = 0.25, 1, and 4 and 
s = 0.3. They will be discussed in conjunction with the non-terminating orbits 
shown in [3j [4] and [5] and with the terminating orbits shown in [6j [7] and [8] in 
Region I when we discuss how the orbital trajectories change as we cross the 
upper boundary of Region I, si, to Region II in Section 6. We shall also discuss 
in Section 6 how the orbital trajectories change as we cross the left boundary 
of Region I, s[, to Region II, and for this purpose, we show an example of a 
terminating orbit in Region II to the left of s[ for e 2 = —0.25 and s = 0.2 in 12 



Thus Region II, with orbits given by eq. (43) , allows only terminating orbits 
with a finite initial distance of the particle from the black hole given by eq.(44) 
for e 2 < 1, and with an infinite initial distance for e 2 > 1. 

We now describe in greater details the boundaries for Regions I and II. 



4 The Boundaries of Regions I and II 

In this section, we present the boundary curves of Regions I and II, and it 
is instructive to see these curves in three different parameter spaces: (e 2 ,s 2 ), 
(k 2 ,s 2 ) and (-g 3l -g 2 )- 

I. The parameter space (e 2 ,s 2 ) or (e 2 ,s) (see [l] and [2| 

There are four relevant boundaries for our universal map in the parameter 
space (e 2 ,s 2 ) for —00 < e 2 < +00: (1) the bottom boundary s 2 = of Region 
I and of Region II (s 2 < is physically inadmissible), (2) the left boundary of 
Region I which we call s — s[, (3) the top boundary of Region I which we call 
s = si, and (4) the top boundary of Region II which we call s — s 2 (we note 
that s > s 2 corresponds to n 2 < and is physically inadmissible). A common 
equation that characterizes boundaries (l)-(3) is A = 0. 

(1) A = 0, s 2 = and k 2 = 0. 
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It is easily verified that setting s 2 — gives g 2 = 1/12, 53 = 1/216, = 0, 
ei = 1/6, e 2 = e 3 = -1/12 from eqs.(14), (22) and (21), and thus A = from 
cq.(17), k 2 = and 7 = 1/2 from eqs.(20) and (19). 

(2) A = 0, k 2 = 0, and s 2 = s' 2 where s' 2 is given by 



_ 1 - 9e 2 - v/(l + 3e 2 ) 3 



* = 27(1 -e 2 ) 2 (48) 

for -1/3 < e 2 < 0, which results in 1/12 > sf > 0. Eq.(48) can be inverted 
to give e 2 in terms of s 2 as 



1 - 18s 2 + 54s 4 - ^(l-12s 2 ) 3 



e , _ ■ - - Ms4 (49) 

for 1/12 > s 2 > 0, which results in -1/3 < e 2 < 0. We shall refer to 
s 2 = s' 2 as the left boundary of Region I. It can be verified that along this 
boundary, 6 = 0, a = (1 - 12s'?) 1 / 2 /6, e 2 = e 3 = -(1 - 12s' 2 ) 1 / 2 / 12, k 2 = 

and 7 = (1 - 12si a ) 1/4 / 2 - 

(3) A = 0, k 2 = 1, and s 2 = s 2 where s 2 is given by 



_ l-9 e 2 + V(l + 3e 2 )3 

' Sl " 27(1 -e 2 ) 2 (50) 

for —1/3 < e 2 < 00, which results in 1/12 > s\ > (with a note that 
s 2 = 1/16 for e = 1). Eq.(50) can be inverted to give 



2 _ 1 - 18s 2 + 54s 4 + ^/(l-^s 2 ) 3 



54s 4 (51) 

for 1/12 > s 2 > 0, which results in —1/3 < e 2 < 00. We shall refer to 
s 2 = s\ as the upper boundary of Region I. It can be verified that along this 
boundary, 6 = it, ei = e 2 = (1 - 12s 2 ) 1 / 2 /12, e 2 = e 3 = -(1 - 12s 2 ) 1 / 2 /6, 
fc 2 = 1 and 7 = (1 - 12«?) 1 /*/2. 

The intersection point of s^ and si is the vertex point (e 2 , s 2 ) = (—1/3, 1/12) 
where the innermost stable circular orbit occurs and where all k 2 = const, curves 
meet. 

(4) s 2 — s\, where s 2 is given by 

4 = ^- (52) 
1 — e z 

We shall refer to s 2 = s\ as the top boundary of Region II. The physical 
requirement that k 2 > 0,where k is given by eq.(6), leads to the condition that 
s 2 < 1/(1 — e 2 ). It can be shown that eq.(52) determines that q 2 , the initial 
position of the particle given by eq.(44), is equal to 1, and thus s 2 > s\ implies 
that the initial position of the particle is inside the Schwarzschild horizon, i.e. 
the physically inadmissible case of q 2 < 1. We call the region s 2 > s 2 Region 
IF. 

II. The Parameter Space (k 2 ,s 2 ) (see 13) 
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This parameter space is equivalent to one used by Hagihara [9, Fig. 13] but 
he did not clearly identify the physical significance of the two parameters. 

There are four relevant boundaries in the parameter space (k 2 , s 2 ) for k 2 > 0, 
s 2 > 0: (1) the bottom boundary s 2 = 0, (2) the left boundary of Region I which 
we call s — s[, (3) the top boundary of Region I which we call s = s 1 , and (4) 
the left boundary of Region II defined by k 2 = 0. 

(1) A = 0, s 2 — and k 2 = 0. This is the lower boundary for Region I and 
for Region II. 

(2) A = 0, k 2 = 0, and s 2 = s' 2 where s' 2 is given by 

, 2 _ -(27k 4 - 36k 2 + 8) - ^k 2 (9k 2 - 8) 3 



32 

for 8/9 < k 2 < 1, which results in 1/12 > sf > 0. It can be inverted to give 
k 2 in terms of s 2 as 



H 2 



1 + 36s 2 - y/(l - 12s 2 ) 3 



54s 4 

for 1/12 > s 2 > 0, which results in 8/9 < k 2 < 1. We shall refer to s 2 = s'i 
as the left boundary of Region I. 

(3) A = 0, k 2 = 1, and s 2 = s 2 where s 2 is given by 

2 _ -(27k 4 - 36k 2 + 8) + v/k 2 (9k 2 - 8) 3 



32 

for 8/9 < k 2 < oo, which results in 1/12 > s\ > 0. It can be inverted to give 



K 2 



1 + 36s 2 + y/(l - 12s 2 ) 3 



54s 4 

for 1/12 > s 2 > 0, which results in 8/9 < k 2 < oo. We shall refer to s 2 = s\ 
as the upper boundary of Region I. 

The intersection point V of s[ and si is the vertex point (k 2 , s 2 ) = (8/9, 1/12). 
Two special k 2 = const, curves, k 2 = 1/2 ± y/3/4 = 0.933012 and 0.0669872, 
are shown as straight horizontal lines to right and left of V (see Section 7, Id). 

(4)k 2 = 0. 

This is the left boundary for Region II that is equivalent to the upper (and 
left) boundary S2 for Region II in the (e 2 ,s 2 ) parameter space. 

Fig. 13 shows these four boundaries in the (k 2 ,s 2 ) parameter space (s 2 = 
and k 2 = are the coordinate axes) that are similar to those in the (e 2 ,s 2 ) 
parameter space. However, from eq.(15), we have 

s 2 = ^-U 2 -l), 
e z — 1 

and thus the constant e 2 curves are straight lines with slopes l/(e 2 — 1) that 
originate from the same point (k 2 , s 2 ) = (1,0), compared to the parallel vertical 
lines in the (e 2 , s 2 ) parameter space. 



III. The Parameter Space (—33,-32) (see 14). 



15 



The reason that we use the coordinates (—33, —32) instead of (32,33) is that 
they can be more easily identified with the parameter spaces (k 2 , s 2 ) and (e 2 , s 2 ) 
so that the horizontal axis —33 increasing to the right corresponds to n 2 (or e 2 ) 
axis increasing to the right, and the vertical axis —32 increasing vertically up 
corresponds to the s 2 axis increasing upward. 

There are four relevant boundaries for our universal map in the parameter 
space (—33, —32) for —00 < —33 < +00, —32 > —1/12: (1) the bottom boundary 
—32 = —1/12, (2) the left boundary of Region I which we call s — s[, (3) the 
right boundary of Region I which we call s = s\, and (4) the left boundary of 
Region II, 33 = gf\ 

(1) A = 0, —32 = —1/12 and k 2 = 0. This is the lower boundary for Region 
I and for Region II. 

(2) A = 0, k 2 = 0, and —33 = —33°^ where 33^ is given by 



(0) '/>\ 3 / 2 



(f) <»> 



for -1/216 < -33 < 0, -1/12 < -32 < 0. We shall refer to -33 = -g { ° } as 
the left boundary of Region I, and is equivalent to s 2 = s' 2 for the (k 2 , s 2 ) or 
(e 2 ,s 2 ) parameter space. 

(3) A = 0, k 2 = 1, and —33 = — where g^ is given by 

^=-(f) 3/2 . (M) 

for < -33 < 1/216, -1/12 < -32 < 0. We shall refer to -33 = -3^ as 
the right boundary of Region I, and is equivalent to s 2 — s\ for the (k 2 , s 2 ) or 
(e 2 ,s 2 ) parameter space. 

The origin (—33, —32) = (0, 0) is the vertex point V where ISCO occurs and 
where all k 2 — const, curves meet. The vertical axis 33 = coincides with 
k 2 = 0.5; the horizontal axis to the right of the origin (—33 > 0) coincides 
with k 2 = 0.933012, and the horizontal axis to the left of the origin (—33 < 0) 
coincides with k 2 = 0.0669872 (see Section 7, Id). 

(4) 33 = 33 2) where gf ] is given by 

(2) 1 1 

This is the left boundary for Region II that is equivalent to k 2 = or 
s 2 = s 2 = 1/(1 — e 2 ) for the (k 2 ,s 2 ) or (e 2 ,s 2 ) parameter space. 

The triangular region bounded by the three curves (1), (2) and (3) given 
above is Region I. The rest of the region above the straight lines given by (1) 
and (4) is Region II. 

While the boundaries in the parameter space appear very symmetrical in the 
parameter space (—33, —32) compared to those in the parameter space (e 2 , s 2 ), 
the e 2 = const, curves are parabolas represented by the equation 
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with the exception that e 2 = 1 is a straight line. These curves all originate 
from the point (-g 3 , -g 2 ) = (-1/216, -1/12). 

The results which are described in the (e 2 ,s 2 ) [or (e 2 ,s)] parameter space 
can be easily translated into the (k 2 , s 2 ) and (—33, —32) parameter spaces using 
eqs.(ll), (14), (15) and the correspondences given in this section. 

5 Orbital Trajectories on the Boundaries 

The orbital trajectories on the boundaries of Regions I and II are important 
special cases. Consider the orbital trajectories given by the three solutions 
given by eqs.(18), (35) and (43) separately. 
Solution (Al) 

The left boundary k 2 = given by s' 2 corresponds to = in eq.(22) and 
e = in eq.(26), and we have e\ = (1 — 12s' 1 2 ) 1 / 2 /6 = — 2e 2 = — 2e 3 . The curve 
for k 2 = and thus s' 2 extends from e 2 = -1/3, s 2 = 1/12 to e 2 = 0, s 2 = 0. 
The orbital trajectories given by eq.(18) on this boundary curve are circular 
orbits the radii q c of which are given by 

or 



Qc = ^ {l + (1 - 12.sf) 1/2 } = ^ {-1 + 3e 2 - (1 + ■Ae 2 ?' 2 } . (58) 

For the values of s' 2 that range from to 1/12 (and -1/3 < e 2 < 0), the 
range of radii for these stable circular orbits is 

3 < q c < 00, (59) 

or 

6GM 

— s— < r c < 00. 
c z 

The circular orbit with the radius q c = 3 (or r c = 6GM/c 2 ) at the vertex V 
represented by the coordinates (e 2 , s 2 ) = (-1/3, 1/12) or (k 2 , s 2 ) = (8/9, 1/12) 
is called the innermost stable circular orbit (ISCO). 

Because every orbital trajectory given by Solution (Al) precesses according 
to eq.(27), these circular orbits also precess with a precession angle 

A0 = 2^{(l-12s' 1 2 )- 1 / 4 -l} (60) 
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for < s'i < 1/12 even though this precession angle is not observable unless 
the orbit is slightly perturbed from its circular path [12]. In terms of the radius 
r c of the circular orbit, the precession angle is given by 

A* = 2*j(l-^) "''-I}. (61) 

Thus excluding the case of zero gravitational field s[ = for which the 
circular orbits have an infinite radius, all (perfectly) circular orbits have a non- 
zero precession angle. The precession angle for the innermost stable circular 
orbit (ISCO), which occurs at the vertex point V is infinite. 

The upper boundary k 2 — 1 given by s\ corresponds to = tt in eq.(22), and 

1 /2 

we have ei = e2 = (l — 12s 2 ) /12 = Substituting these into eqs.(23), 

1/2 

(24) and (26) gives a useful relationship (l — 12s 2 ) = 2e/(3 + e). The curve 
for k 2 = 1 given by s\ extends from e 2 = —1/3, s 2 = 1/12 to e 2 = oo, s 2 = 0. 
Equation (18) on this boundary becomes 

i- 1 | ( 1^5secM2^) \ 

9 _ 3 + 3V 1 12s i 1, l + S ec/i(2 7 0) J' ( J 

where 

7 =^(1-12 S 2 ) 1/4 . (63) 

For —1/3 < e 2 < 1 (1/12 > s 2 > 1/16), the particle starts from an initial 
position q max at 4> — given by 

and ends up at = oo circling the star or black hole with a radius that 
asymptotically approaches g m ; n given by 

For e 2 = 1, the particle starts at infinity at the polar angle <f> = 0, and for 
e 2 > 1 (s 2 < 1/16), the particle starts from infinity at the polar angle <j> — 
given from eq.(32) by 

tanh 2 ( 7 * 1 ) = ^-^(l-12 S ?)- 1/2 , (66) 

and ends up at <f) — 00 circling the star or black hole with a radius that 
asymptotically approaches q min given by eq.(63). The angle can range from 
for e 2 = 1 (parabolic-type orbit) to 2tanh~ 1 (l/y3) = 1.31696 = 75.456° for 
very large e 2 . 

We refer to all these orbits given by eq.(62) as the asymptotic orbits. 
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Examples of asymptotic orbits for e 2 = 0.25, 1 and 4 are shown in 15 16 
and El 

Solution (A2) 

Substituting k 2 = and s 2 — s' 2 into eq.(35) gives the equation for the 
orbital trajectories for Solution (A2) on the left boundary of Region I as 



where 



1 _ 1 

q ~ 3 



1 - (1 - 124' 



(1 - 12si 2 
cos 2 (7^) 



1/2 



(67) 



1 



7 



(1 



i2 S ; 2 ) 1/4 . 



(68) 



Thus on the left boundary of Region I, the terminating orbit (inside Region I) 
given by Solution (A2) remains a terminating orbit. The particle, starting from 

1/2 

the polar angle </> = at a distance qi given by 1/q-y = 1/3 + 2 (l — 12s' 2 ) /3 

plunges into the center of the black hole when <j> = <j>i = n/ (l — 12s'] 2 ) . 

On the other hand, substituting k 2 — 1 and s 2 = s\ into eq.(35) leads to q = 
const, independent of 4> (because e\ = e-i and the dependence on <j> cancels out). 
Thus on the upper boundary of Region I the terminating orbit (inside Region 
I) given by solution (A2) becomes an unstable circular orbit of radius q u given 
by 



(1 - 12s 2 



2\l/2 



2{l-3e 2 + (l + 3e 2 ) 1 /2} _ 1 + £ 



9(1 



3 + e' 



(69) 



and 



if {l - (1 - 12*?) 172 } = ^2 {-1 + 3e 2 + (1 + Ze 2 ) 112 } , (70) 



4*i 



Qu 



- t -r <71) 

It is an unstable circle because its radius q u is identical to the radius <7 m i n of 
the asymptotic circle given by eq.(65). For the values of s 2 that have a range 
< s 2 < 1/12 ( and oo > e 2 > —1/3), we have the range of values for the radii 
of these unstable circular orbits along the upper boundary s 2 — s 2 of Region I 
given by 

1.5 < q u < 3, (72) 



2>GM 6GM 
— s— <r u < — ^— ■ 
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Notice the similarity between the expressions given by eqs.(69) and (70) for 
the radii of the unstable circles on si and those for the radii of the stable circles 
on s[ given by eqs.(57) and (58). Unlike the stable circles for which e = 0, these 
unstable circular orbits do not have e = except when q u — > 3 which is the 
radius of the innermost stable circular orbit (ISCO) at the vertex V. 

Solution (B) 

Substituting k 2 = and s 2 = s' 2 into eq.(43) gives the equation for the ter- 
minating orbits for Solution (B) on the left boundary of Region I. This equation 
coincides with eq.(67) which is the expression that describes the terminating 
orbits for Solution (A2) on the same boundary i.e. with q 2 of eq.(44) becoming 
qi of eq.(36) and with a = e\ = (1 — 12s' 1 2 )/6 on the boundary. 

Substituting k 2 = 1 and s 2 = s\ into eq.(43), the terminating orbits for 
Solution (B) on the lower boundary of Region II (which is the upper boundary 
of Region I) become the same asymptotic orbits given by eq.(62) which describes 
the non-terminating orbits of Solution (Al) on the same boundary. That is, on 
the boundary s 2 = s\ , q 2 of eq.(44) becomes <7 max of eq.(64) and with a — e 3 = 
-(1 - 12sf)/6. The value of q mm is given by eq.(65) for —1/3 < e 2 < 1. When 
e 2 > 1 the particle starts from infinity (q 2 = oo) at a polar angle <fi = that 
is given by eq.(66). 

Thus for the orbits on the boundaries given by Solutions (Al), (A2) and 
(B), we have (1) on s[, the stable circular orbits of radii q c given by eq.(58) 
and the terminating orbits given by eq.(67), and (2) on si, the unstable circular 
orbits of radii q u given by eq.(70) and the asymptotic orbits given by eq.(62). 
Even though these results reproduce the conclusions obtained from the effective 
potential calculation, they place the orbits and their corresponding analytic 
expressions clearly on the boundary curves and si shown in [2] 

6 Trajectory Changes Across the Boundaries 

. 2 

The left hand side of eq.(7) excluding the r term can be considered to be an 
effective potential. In a similar manner we define, from eq.(13), an "effective 
potential" W to be 

W = ~(4U 3 - g 2 U - g 3 ). 

Although only the condition that W < is physically realizable, it is useful 
to see the entire curve W above and below the horizontal U axis and we shall 
examine it and its variation in this section. 

We want to show the changes in the orbital trajectory in conjunction with 
the changes of this effective potential as we consider various coordinate points 
on a counterclockwise path in parameter space around the vertex point V in 
[2j We have chosen a specific path that generally represents the changes in the 
trajectories and in the effective potential as boundaries are crossed. The points 
where the W curve intersects the U axis are the real roots of eq.(16). Recall 
that Solution (Al) for Region I requires e\ > e 2 > U > es, Solution (A2) for 
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Region I requires U > e\ > e 2 > e 3 , and Solution (B) for Region II requires 
U > a. For U denned by eq.(12), U — > +00 implies r — > 0. 

We start by considering a coordinate point inside Region I for which e 2 < 1 
that we choose to be (e 2 ,s) = (0.25,0.2) (see [2]). We know from Section 3 that 
it allows a precessional elliptic- type orbit given by eq.(18) shown in [3j and a 
terminating orbit given by eq.(35) shown in Fig. 6. Note that q\ is less than 
q m in and the two possible orbits describe two separate initial conditions. The 
corresponding W curve is shown in 18 Noting that e 3 > — 1/12 for this example, 
the bowl-like "potential" below the U axis between U = e 3 and U = e 2 allows 
an elliptic-type orbit for a particle between q = q max and q = q m [ n given by eqs. 
(23) and (24). On the other hand, the down-hill "potential" below the U axis 
for U > e\ allows a terminating orbit for a particle that starts at U — e\ and 
ends at U = +00 or from q = q\ given by eq.(36) to q = 0. Thus the possible 
orbits shown in [3] and [6] are appropriately matched with the W curve. 

Going counterclockwise around the vertex V (see [2]), we choose the next 
coordinate point (0.25,0.26481) which is straight above the previous point but 
is on the boundary S\ characterized by k = 1. This point allows an asymptotic 
orbit given by eq.(62) and shown in 15 with g max and <7 m i n given by eqs. (64) 
and (65), and it also allows a circular orbit with a radius q u given by eq.(69) 
with q u = q m i n . The corresponding W curve is shown in [19] which indicates 
that ei = e 2 . The bowl- like potential between U = e 3 and U — e\ = e 2 
allows a particle starting from U = 63 or q = g max to circle the black hole 
with a radius that asymptotically approaches U — ei = e 2 or q = q nim . The 
down-hill potential for U > e\ — e 2 allows the particle to remain (unstably) at 
U = e\ = 62 or for it to have a circular orbit of radius q = q u . This circular 
orbit is clearly unstable because a small perturbation would make it into a 
terminating or an asymptotic orbit. 

We next choose a point (0.25, 0.3) in Region II which is again straight above 
the previous point. It allows only a terminating orbit given by eq.(43) and 
shown in |9j The corresponding W curve is shown in [20] with one real root at 
U — a. The transition from Region I to Region II across Si clearly shows that 
e3 becomes a. Noting that a > —1/12, the down-hill potential for U > a allows 
a particle that starts at a finite distance from the black hole at U — a or q = q 2 
given by eq.(44) and ends at U — +00 or q = 0. The asymptotic orbit shown in 
|15| becomes a terminating orbit shown in [9] as we go from the boundary s± to 
Region II. The unstable circular orbit on si also becomes a terminating orbit (as 
the initial distance becomes less than q u ) which is part of the same terminating 
orbit in Region II arising from the asymptotic orbit on the boundary (see the 



transition of the W curve below the horizontal axis from 19 to 20). These 



examples show how Solutions (Al) and (A2) merge into Solution (B) across the 
boundary s±. 

We choose the next coordinate point (—0.5,0.3) in Region II that is hori- 
zontally left to the previous point (see [2]). The W curve is shown in 21 and 
indicates the one real root U — a has increased in the positive direction, i.e. the 
initial distance q 2 of the particle from the black hole has decreased compared to 
the previous point. The terminating orbit is similar to [9] except that the initial 
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distance is smaller. If the chosen coordinate point were to be further left in[2]so 
that it fell on the left upper boundary S2 of Region II (see [I]), a would be equal 
to 1/6 and the initial position of the particle would be on the Schwarzschild 
horizon corresponding to q<x = 1- 

We next select the coordinate point (—0.25, 0.2) in Region II below the previ- 
ous point but close to the left boundary s[ of Region I (see [2]). The terminating 
orbit is shown in Fig. 12. The W curve is shown in Fig. 22 where the real root 
has become more positive while the left part of the curve is slightly above and 
almost touches the U axis. If the chosen coordinate point were to be further 
down in Fig. 2 so that it fell on the lower boundary s — of Region II (see Figs. 
1 and 2), a would be equal to 1/6 and the initial position of the particle would 
be on the Schwarzschild horizon corresponding to q% = 1. 

We now go to the right and choose the next coordinate point (e 2 ,s) = 
(—0.099275,0.2) exactly on the boundary characterized by k 2 — (see [2]). 
The W curve is shown in [23] It touches the U axis at U — e-i — e 3 and allows a 
circular orbit with a radius q c given by eq.(57). The down-hill potential given by 
U > e\ gives a terminating orbit with q\ given by eq.(36). The circular orbit is 



a stable one as a small perturbation that could change [23\ into one similar to 18 



clearly confines the perturbation to be bounded, and the possible terminating 
orbit solution is well separated from it as e\ is greater than = 63, or qi is less 
than q c . The terminating orbit for this point is similar to the terminating orbit 
of the previous point except for a small difference in q± . The transition of the 
terminating orbit from Region II to Region I is a smooth and gradual one (with 
a becoming e{). However, when we cross the boundary s[, a stable circular orbit 
becomes possible on the boundary and it then becomes an elliptical orbit as the 
coordinate point moves into the interior of Region I. As we noted previously 
[see eqs. (60) and (61)], the circular orbits, like the elliptical orbits, precess. 

We now return to our starting point (0.25, 0.2) for which the W curve is 
shown in 18 and the possible trajectories are shown by [3] and [6] 

In the plots for the "potential" curve W, it is useful to mark the U axis 
where U < —1/12 by a thick line for the following identification. As we have 
already noted, Solution (Al) gives orbits that are elliptic-type, parabolic-type 
and hyperbolic-type for e% >,=, and < —1/12 respectively, and Solution (B) 
gives terminating orbits for a particle that is initially at a distance q^ from the 
black hole that is finite, infinite at <fi = 0, and infinite at ^2 > for a >,= 
and < —1/12 respectively. The examples which illustrate the crossing of the 
si boundary that are shown by [3} [T5] and [9] for the orbital trajectories, and by 



18 19 and 20 for the W curves, are for the cases and a greater than —1/12 
that are representative of e 2 < 1. If we chose to cross the Si boundary along a 
path on e 2 = 1, the W curve would have = —1/12 or a = —1/12, and the 
corresponding representative trajectories would be like those shown in Figs. 4, 
16 and 10; and if we chose to cross the si boundary along the path on e 2 > 1, 
the W curve would have < —1/12 or a < —1/12 (i.e. the left side of the W 
curve would cut the thick line representing U < —1/12) and the corresponding 
representative trajectories would be like those shown in [5j [17] and 11 As for 



the terminating orbits in Region I, they are not affected much by the values of 
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e 2 , as are shown by[6j [7] and [8] and when we cross the boundary s\, they all 
become unstable circular orbits on the boundary and then become parts of the 
terminating orbits in Region II. 

We see that the shape of the W curve and whether or not it crosses the 
thick line that represents U < —1/12 can be used to foresee the possible type 
of orbits and trajectories and to relate them to the approximate location of the 
coordinate point on the map. 

The W curve at the vertex V located at (e 2 , s) — (—1/3, 1/VT2) is given by 
W = — 4t/ 3 and it has three real roots coinciding at e\ — e 2 = e 3 = 0. 

7 Analytic Expressions for Special Cases 

The exact analytic orbit equations (18), (35) and (43), and their corresponding 
equations on the boundaries, eqs.(57), (62), (67) and (69), give us a clear picture 
of what these orbits look like and what their coordinates are on a map char- 
acterized by (e 2 ,s 2 ) and with boundaries of Regions I and II characterized by 
eqs.(48), (50), and (52). While generally we are required to start the numerical 
computation of the orbits with eq.(21) or (39), there are a number of special 
values of squared modulus k 2 and squared energy parameter e 2 that allow var- 
ious exact relations that are very interesting and useful. We have exploited the 
theory of Jacobian elliptic functions to come up with many of these relation- 
ships. A few of these analytic relations have been presented in a previous paper 
[2], but it is convenient to collect all of them in one place and we present and 
summarize them in the following. 

(I) Curves of k 2 = const, and the Orbit Equations on These Curves. 

(la) For k 2 = 0, Sections 4 and 5 give many special results for this special 
k 2 value. In particular, we have eq.(48) [or (49)] that gives the coordinates 
(e 2 , s' 2 ) of this curve. The radii of the stable circular orbits from Solution (Al) 
are given by eqs.(57)-(59) and the precession angle is given by eq.(60) or (61). 
The terminating orbit equation from Solution (A2) is given by eq.(67). 

(lb) For k 2 — 1/2, we have a simple analytic expression that relates the 
values of s 2 and e 2 on it: 



for e 2 > —1/2, < s 2 < 1/9 (in Region I and part of Region II) and for 
which s 2 = 1/18 for e 2 = 1; and a simple expression 



for 1 > e 2 > -1/2, s 2 > 1/9 (in Region II). Equations (73) and (74) can be 
inverted to give 




(73) 




(74) 
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Equations (73) and (74) or eq.(75) are shown as a single continuous curve 
marked by k 2 = 0.5 in Fig. 2. In the parameter space (k 2 , s 2 ), the k 2 = 0.5 curve 
is given by 

S 2 = 1 



18(3k 2 - 2)' 

In the parameter space (—33, —52), the k 2 = 0.5 curve is simply the vertical 
axis 53 = 0. 

By noting that in Region I, e\ = — e 3 = e 2 = 0, 7 = tfgi, the orbit 

equation (18) becomes 

- = \-2^- 2 cn 2 (^,l/V2), (76) 
q 6 

where the values of s 2 in the range 1/12 > s 2 > for g 2 are given by eq.(73), 
and q m i n = 3 independent of e 2 . For the elliptic-type orbits corresponding to 
— 1/3 < e 2 < 1, the initial distance <7 max at (f> = of the particle from the star 
or black hole is given by 

-^^-2732. (77) 

For the parabolic- type orbit corresponding to e 2 = 1, g max = 00 at <p = 0. 
For the hyperbolic- type orbit corresponding to e 2 > 1, the particle approaches 
from infinity along an incoming asymptote at an angle given by eq.(29) to 
the horizontal axis, goes around the star or black hole, and leaves along an 
outgoing asymptote at an angle where 



*' = ^"~V~^W (78 ' 

and 

9i _**em_ 9u (79) 



where K(l/y/2) = 1.85407. 

The terminating orbit equation (35) in Region I given by Solution (A2) 
becomes 



1 1 2 ^ (80) 



q 3 cn 2 ( 7 0, 1/^2)' 
The particle, starting from the initial distance q\ given by 
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- = \+^V9~2 (81) 
plunges into the center of the black hole when its polar angle <j>\ is given by 

K(l/V2) 



n 



K.I2 



(82) 



For the terminating equation (43) in Region II given by Solution (B), the 
values of s 2 are given by eq.(73) for the range —1/3 < e 2 < —1/2, 1/12 < s 2 < 
1/9, and are given by eq.(74) for the range e 2 > —1/2, s 2 > 1/9. By noting from 
eqs.(39)-(42) that A = —B = Vl2s 2 - 1/12, a = 0, 7 = [(12s 2 - l)/3] 1/4 /2, the 
orbit equation becomes 

( 1 ~ <-ii\2-<.>. :i/\'2) 1 ( . s t 



9 3 V 1 - cn(270, 1 



The A; 2 = 1/2 curve in Region II is confined to e 2 < 1 
1/12 to 00, and thus the initial distance of the particle from the black hole is 
always bounded. Indeed qi = 3 independent of e 2 . The particle, starting from 
this distance from the black hole, plunges into the center of the black hole when 
its polar angle </>2 is given by 

7 

The above results supplement the results given in Appendix C of ref. [2] that 
have a gap in the range —1/3 < e 2 < 0. 

(lc) For k 2 — 1, Sections 4 and 5 give many special results for this special k 2 
value. In particular, we have eq.(50) [or (51)] that gives the coordinates (e 2 , s 2 ) 
of this curve. The non-terminating orbit equation given by Solution (Al) 
becomes an asymptotic orbit eq.(62). The relevant quantities for —1/3 < e 2 < 1 
are given by eqs.(64) and (65), and for e 2 > 1 are given by eq.(66) and the 
description after it. The terminating orbit given by Solution (A2) becomes an 
unstable circular orbit with the radius given by eq.(69) or (70) that has a range 
given by eq.(72). The terminating orbit in Region II given by Solution (B) 
becomes an asymptotic orbit given by eq.(62). 

(Id) k 2 = l/2±v / 3/4 are special values because they correspond to s 2 = 1/12 
or <?2 = 0. Substituting s 2 = 1/12 into gs and k 2 given by eqs.(14) and (42) 
show the following: k 2 = 1/2 + \/3/4 = 0.933012 for all values of 33 > which 
holds when k 2 < 8/9 or when e 2 < -1/3; and fc 2 = 1/2 - V^/4 = 0.0669872 
for all values of g% < which holds when k 2 > 8/9 or when e 2 > —1/3. Thus 
the k 2 = 1/2 ± V3/4 curves are represented by a single horizontal straight line 
(s 2 = 1/12 or 32 = 0) and passing through the vertex point V in all our three 



parameter spaces shown in [2| [13] and 14 



(II) Geometric Eccentricity e Defined by Eq.(25) for —1/3 < e < 1 in 
Region I. 
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(2a) £ — > e as s — > for < e < 1; and e = 1 when e = 1, i.e. e = e = 1 for 
all values of s (see ref.[2]). 
(2b) e = for k 2 = 0. 

(2c) For k 2 = 1/2, the following equations relate the values of s 2 and e 2 on 
the curve given by eq.(73) to the value of e defined by eq.(25) for non-terminating 
orbits in Region I in the range — 1/3 < e 2 < 1. 

2 _ i 12(l + £) 2 (l-£)(l + 3 £ ) 

6 - 1 (3 + 6e - e 2 ) 2 ' (85) 

or 



-3(1 + e 2 ) + 2^3(1 + 2e 2 ) + 2\/l - e 2 ^-3(l + e 2 ) + 2^3(1 + 2e 2 ) 

£_ 7-e 2 -2 v / 3(l + 2e 2 ) ' 

(86) 

with e = 1 when e = 1; and 

3 + 6e- £ 2 

~ 36(1 + e) 2 ' 180 

or 

_ 3(1 - 12s 2 ) + 2^/3(1 - 12s 2 ) 
£ _ 1 + 36s 2 ' 

(2d) For k 2 — 1, the following equations relate the values of s 2 and e 2 on 
the curve given by eq.(50) or (51) to the value of e defined by eq.(25) for non- 
terminating orbits in Region I in the range — 1/3 < e 2 < 1. 



e 2 = ( 1+£ ^ 3 + 5£ ) (89) 
(3-e) 2 



or 



and 



or 



e= -(l + 3e») + Vr + ^ 



_ (3- £ )(l + g ) 
" 4(3 + £) 2 ' 1 j 



1 - 12s 2 + 2-y/l - 12s 2 



(92) 



1 + 4s 2 

(III) The Special Case of e 2 = 1. 

(3a) For < s 2 < 1/16 in Region I. The k 2 and 7 given by eqs.(20) and (19) 
for the non-terminating and terminating orbit equations (18) and (35) can be 
expressed simply in terms of s 2 by the following equations: 
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or 



, 2 _ 1 - 8s 2 - VI - 16s 2 
1 2s 2 



(94) 



7 ' 4(1 + fc 2 ) l- Vl-16s 2 ' (95) 
For the non-terminating orbits, q max = 00, and 

1 + fc 2 fQR . 

<7min = fc2 ■ (96) 

For the terminating orbits, the initial distance qi of the particle from the 
black hole is 

qi = l + k 2 , (97) 

and the particle plunges into the center of the black hole, for < k 2 < 1, 
when its polar angle <j> is equal to 4>\ given by eq.(37). For the particular case 
of k 2 — 1, the terminating orbit becomes a circle of radius q c = q\ = 2. 

(3b) For s 2 > 1/16 in Region II. The k 2 and 7 given by eqs.(42) and (41) 
have the following simple expressions: 

7 2 = ^ (99) 

The particle, coming from infinity at <f> = 0, plunges into the center of the 
black hole at a polar angle <fi — <t>2 given by eq.(45). 

(IV) The Special Case of e 2 = -1/3. The 7, k 2 and q 2 given by eqs.(41), (42) 
and (44) for the terminating orbits in Region II have the following expressions: 

7 2 - -±= {(1 - 12,s 2 ) 2 / 3 [l + (1 - 12s 2 ) 1 / 3 + (1 - 12s 2 ) 2 / 3 ]} 1/2 , (100) 



2 _1 V3 (l-12s 2 )V 3 l + (l-12s 2 )V 3 

2 4 1(1 - 12s 2 )V3| [l + (l-12s 2 ) 1 /3 + (l_12 S 2)2/3]l/2' ^ > 



and 



1 1 

?2 



I {l + (1 12s 2 ) 1 / 3 + (1 - 12s 2 ) 2 / 3 } . (102) 
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These expressions for j 2 , k 2 and q 2 follow by substituting into eqs. (39)-(42) 
and (44) A = s 4 (l - 12s 2 ) 2 /(2 4 • 3) for e 2 = -1/3. The particle, initially at a 
distance q 2 from the black hole, terminates its orbit at the center of the black 
hole at a polar angle = <j) 2 given by eq.(45). 

(V) The Maximum Value of q 2 In Region II For A Given Value of e 2 < 0. 

For any given value of e 2 < —1/3, as we increase s 2 from to oo, q 2 in 
Region II given by eq. (44) increases from 1 at s 2 = to a maximum value given 
by 

1 - e 2 

92 max = 5- (103) 

which occurs at a value of s 2 given by 



s 2 = 



-e 



2 



2(1 -e 2 ) 2 ' (1 ° 4) 

and as s 2 increases beyond this value, q 2 decreases and becomes equal to 1 
when s 2 = s\ given by eq.(52), and q 2 continues to decrease to as s 2 increases to 
00. The same expressions for q 2 max and the value of s 2 where q 2ma x occurs apply 
for —1/3 < e 2 < except that we need to exclude the values of s' 2 < s 2 < s\ 
in Region I for which q 2 is not defined. On the boundary of Region I and of II, 
the value of q 2 is given by l/q 2 = 1/3 + 4a, where a = (1 — 12s' 1 2 )/6 on s' 2 and 
a = -(l- 12s 2 )/6 on s 2 . 

Equations (103) and (104) can be proved as follows. We note that q 2 is 
given by eq.(44) and that a satisfies eq.(16), i.e. 4a 3 — g 2 a — — from which 
setting da/d(s 2 ) — gives a = —1/12 + (1 — e 2 )s 2 /2. Substituting this value 
of a into the cubic equation satisfied by a gives the expression for s 2 given by 
eq.(104) which in turn gives a simple expression for a and the simple expression 
for 92 max given by eq.(103). 

The 7 2 and k 2 for the s 2 given by eq.(104) where q 2 = q 2 max can be obtained 
from eqs. (41) and (42) and have the following expressions: 



^ 2 = w^-y (105) 



and 



.2 1 . l + 2e 2 



fc -2 + 4(^)' ( 106 ) 

from which the polar angle <p 2 with which the particle plunges into the center 
of the black hole can be calculated from eq.(45). 

While the descriptions of the orbital trajectories in Regions I and II already 
provide a clear picture of what a particle's orbit would look like in different 
parts of the map (see [l] and |2|, the various analytic formulas given in this 
section can be used to provide quantitative checks of the parameters associated 
with the orbits at various specific coordinate points of the map. We mention 
three examples: (i) the values of k 2 given by eq.(94) along e 2 = 1 in Region I 
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(from s 2 = to 1/16), and those given by eq.(98) in Region II (from s 2 = 1/16 
to oo) provide useful checks on how the various k 2 — const, curves radiating 
from the special vertex point V cross the e 2 = 1 line, (ii) as can be checked 
from eqs.(103) and (104), the values of of eq.(44) increase from 1 to 4 along 
e 2 = —1/3 as s 2 increases from to 3/32, and decrease from 4 to 1 as s 2 
increases from 3/32 to 3/4 (the upper boundary of Regions II), and (hi) the 
geometrical eccentricity e given by eqs.(85)-(91) for the precessing elliptic-type 
orbits corresponding to the specific k 2 = 1/2 can be easily tracked as it varies 
from e = 1 for e 2 = 1, s 2 = 1/18 to e = (2/^-3 + 2^3 - l)" 1 = 0.516588 for 
e 2 = 0, s 2 = (3 - \/3)/18 = 0.0704416. 




8 Newtonian and Non-Newtonian Correspondences 

We have used e 2 defined by eq.(15) as one of the two principal parameters for 
representing the orbit. In this section, we discuss the relation of this parameter 
to the corresponding Newtonian parameter. It will be seen that some relativistic 
bounded orbits do not have corresponding Newtonian orbits when appropriate 
limits are taken. 

From eqs.(15) and (7), the squared energy parameter e 2 can be expressed as 

<?=\^\ {r + \r*-— \ -^ 2 |. (107) 

We first show how in the Newtonian limit the non-terminating orbits in 
Region I given by eq.(18) become the Newtonian elliptical, parabolic or hyper- 
bolic orbit and how the energy parameter e becomes the "eccentricity" of the 
Newtonian orbit. 

In Newtonian mechanics, the energy E of a particle of mass m in a gravi- 
tational field produced by a massive body of mass M is the sum of its kinetic 
and potential energies and is given by 

E = 1 -m{r\r^)-^, (108) 

where the derivative ■ represents d/dt, t being the ordinary time. If we define 
the Newtonian eccentricity ex to be 



1/2 

ZN = \ 1 + ^*V2\ > ( 109 ) 



2E h 2 
m{GMf 



. 2 

the Newtonian equation of motion r—rcp = -GM/r 2 obtained from eq.( 108) 
can be integrated to give the equation for the orbit given by 

1 = ^(l-e w cos0). (110) 
r 
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Equation (110) is used with the usual identification that < ejv < 1 gives 
an elliptic orbit, e^v = 1 gives a parabolic orbit, and cn > 1 gives a hyperbolic 
orbit. Note that e^v is defined for the range of all positive values < ejv < oo. 

Before we discuss how eqs.(107) and (109) are related, and how eq.(llO) is an 
approximation of eq.(18), it is useful to review briefly the curves of conic sections 
purely from the mathematical point of view. It is known that the equation in 
the polar coordinates (r, <j>) 



represents a curve traced by a point, say P, about a focus, say F, on the 
x-axis, such that the ratio of the distance r = PF to the perpendicular distance 
PD to a directrix which is a line parallel to the y-axis, is a constant equal to 
£c, i-c sc = PF/PD, and ec is defined for all positive values and thus in 
the range < ec < oo. The parameters ec and p characterize three types of 
curves. The curve is an ellipse for < ec < 1 and p = a(l — e c ), a parabola 
for ec = 1 and p = 2r m ; n , and a hyperbola for ec > 1 and p — a(e c — 1). For 
the elliptic curve, the parameters a and b = a^J\ — e c are associated with the 
semi-major and semi-minor axes, and ec can be shown to be related to r m - m 
and r max by e c = (r max - r min )/(r max + r min ). This expression for e c in terms 
of r max and r m j n is not valid for the hyperbolic curve. For the hyperbolic curve, 
the parameters a and b — ay/ ' e c — 1 are used to characterize the angles with 
respect to the x-axis of the two asymptotes. 

By identifying the equation of curves of the conic sections given by eq.(lll) 
with the Newtonian orbit equation (110), we say that the Newtonian orbit is 
an elliptic, parabolic or hyperbolic curve depending on the value of e^v, and we 
identify e» with ec for the entire range of positive values from to oo, and by 
equating GM/h 2 = 1/p. 

We now want to see what the Newtonian limit implies from our general 
relativistic orbit equation (18). First, from eqs.(108) and (109), the Newtonian 
eccentricity ejv can be expressed as 



since setting c = oo is not consistent with reality, we will proceed by stating 
that 6 — y 6 7v if we take the approximation r — > t and 



Next, we proceed from eq.(18) to eq.(llO) in steps which, by noting eqs.(107) 
and (113), must be taken with care around the point e 2 = 0. For small s 2 , 



- = -(1 — eccoscb) 
r p 



(111) 





(113) 
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expanding the parameters and functions in eqs.(18)-(22) in power series in s, 
we find [3], for e 2 > 0, 



ei = - - s 2 - (e 2 + 3)s 4 + .. 

1 !/ -x 2 (e + 1) 3 4 
62 = "12 + 2 (e + 1)5 + + 



e3 



1 1 

12 ~ 2 



(e - l)s 2 + 



2e 

(e-1) 3 
2e 



(114) 



A: 2 = 4es 2 1 + 



l + 9e 2 -2e 3 2 



7 =-[l-(3-e)s 2 



(115) 
(116) 



if(fc) = -(l + es 2 + ..), 



(117) 



and 



sn(jcf), k) = sin - 



e sin 



sin ■ 



3<A 3<j> 



„ cos — 
2 2 



a' + ... 



Substituting the above approximations into eqs.(18), (26) and (27) gives, to 
the order of s 2 , 



= 2s 2 {l-ecos[(l-<5)0]}, 



which, in terms of r, is the approximate orbit equation 

£ = ^{l-ecos[(l-5)fl}, 
where e, to the order of s 2 , is given by 



e ~ e + (e^ 1 



e 3 )s 2 , 



and where 5, to the order of s 2 , is given by 

6 ~ 3s 2 . 



(118) 



(119) 



(120) 



(5 is related to the precessional angle given in eq.(27) by A<j> ~ 2-7r(5 ~ 
67rs 2 = 67r[GM/ (/ic)] 2 and it is independent of e (to the order s 2 ) for < e < oo. 
Thus if we ignore terms of order s 2 , we recover the Newtonian orbit equation 
given by 



GM 



(1 — ecos^). 



(121) 
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The question that can be posed at this point is whether we can define the 
Newtonian limit by stating that it is the general relativistic result for small s if 
we ignore terms of order s 2 and higher. 

To proceed, we note that the case e = is excluded from the expansions 
given by eqs.(114)-(117) and it is also clear that e 2 = does not satisfy the 
condition given by eq.(113) and must be treated separately. 

For e 2 = 0, instead of eqs.(114)-(117), we find 



ei 



s 2 - 3s 4 - 20s 6 + ... 



e2 



e 3 = 



1 2 , 3 3 



12 



s b + 10s h + ... 



1 1 9 v2 o 3 d 21y2 r 
~12 + 2 S + 2 S ~^T S +1 ° S + "" 



(122) 



4V2 S 3 [l + ^s 2 + ...), (123) 



1 3 2 V2 3 27 4 33^2 5 
7=2-^ + "T* + — S +- (124) 

K{k) = |(l + \/2s 3 + ...) , (125) 



and 



sn(70, fc) = sin - - ys cos - + — ^sin - + sin — j s + ... 

Substituting the above approximations into eqs.(18), (26) and (27) gives the 
orbit equation for e 2 = in the form given by eq.(118) with the same 5 given 
by eq.(120) but, instead of cq.(119), e is now given by 

s = V2s (l + ^s 2 ^j (126) 

to the order of s 2 . Thus we have an elliptic orbit that precesses with the 
same angle S given by eq.(120) but with an eccentricity equal to a/2s + 9\/2s 3 /4 
to the order s . If we ignore terms of order s 2 , the orbit equation becomes 

1 GM , r , x 

- = -p-(l-V^*coB^), (127) 

which is a (non-precessing) elliptical orbit with a non-Newtonian eccentricity 
that is dependent on the speed of light. To the best of our knowledge, this 
general relativistic elliptical orbit with an eccentricity of the order s has not 
been noted by other authors. As an example, an elliptic orbit with eccentricity 
£ = 0.079005 could have resulted from the general relativistic orbit with e = 
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[remembering that e is denned by eq.(107) and not eq.(112)] and s = 0.0554787 
for which the approximation formula (126) gives e ~ 0.079002 and the first 
term alone gives s/2s = 0.078459, and eq.(120) gives 5 ~ 0.009234. The orbit 
for e = from general relativity becomes a Newtonian circular orbit if we ignore 
the first order term in s, i.e. ignoring the second order term s 2 and higher order 
terms is not sufficient to get a Newtonian limit for this case. 
It is clear that the entire region characterized by e 2 < or 

.2 / ■ GmV 2GM 2 ... 
r + [r<j> < — (128) 

V r 2 cj>J c 

is non-Newtonian in character. This includes all bound orbits in Region I 
in the range —1/3 < e 2 < 0, and all circular orbits that occur on the curve s' 2 
from the vertex V to the origin O for which k 2 — and s = 0. All circular 
orbits precess even though the precession angle is not observable [12], and the 
precession angle is given by eq.(61) which is non-zero unless the radius of the 
circle is infinite which occurs on = for zero gravitational field. For small s, 
the precession angle is given approximately by 

. j „ 2 GM 
A(p ~ bns ~§ir — — , 
c z r c 

where r c ~ h 2 /{GM) is the radius of the circular orbit. For small sf, the 
values of e 2 along the s[ curve where the circular orbits occur are given by 

e 2 ~ -2sf. 

For small and for s just above inside Region I given by s 2 = s' 2 + (As) 2 , 
it can be shown that to the order As we have an elliptic orbit similar to eq.(127) 
given by 

" = ^[l-\/2(As)cos</>]. (129) 
r h z 

To summarize the above results, for small values of s, the general relativistic 
correction to the Newtonian elliptic orbit is second order in s for e 2 > for 
which the correction appears in the angle of precession, but is first order in s or 
As for e 2 < for which the correction appears in the eccentricity of the orbit. 
This division gives a new meaning and significance to the parameter e 2 . The 
new feature of this result is that if we ignore terms of order s 2 and higher but not 
terms of order s, then there exist non-Newtonian elliptic orbits of eccentricity 
given by eq.(127). 

For the hyperbolic-type orbit (e 2 > 1), in the Newtonian limit for small s 2 , 
eq.(28) becomes 

sin 2 ^^ 1 --^, (130) 
2 2 2e 1 

from which we find 
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l-i-cos- 1 ^)^, (131) 

and the complementary angle = 2ir — ^2, where ^2 is given by eq.(30), 
becomes ^ ~ 0o also. The Newtonian hyperbolic equation is thus given by 
cq.(llO) or (121) with e N = e > 1 where cos" 1 ^/ 6 ) < < 2tt - cos-^l/e). If 
we define 

Q GR = *! + *i (132) 

and 

©TVetoton = 200 , (133) 

the difference A0 = &Newton — Qgr can be taken to be an analog of the 
precession angle given by eq.(27) for a hyperbolic orbit, and for small s 2 , it was 
shown [3] to be given by 

A0 ~ [6tt - 60o + 2(2 + e" 2 ) Ve 2 - l] s 2 . (134) 

This result is different from the approximation used by Longuski et al. [13]. 
As discussed in ref.13, an experimental test can be carried out to check this 
result. 

We see that our parameter e 2 defined by eq.(15) or (107) not only has the 
advantage of being directly associated with the total energy but is also directly 
related with e 2 N defined by eq.(109) or (112) in the Newtonian limit for the 
range < e 2 < oo. The region characterized by e 2 < alerted us that there is 
a region other than the high gravitational field (s 2 >> 0) region in which the 
orbits are non-Newtonian in character; these include all those orbits for which 
e 2 < 0, and all stable circular orbits along for all values of s. In particular, 
for small s and to the order s, a non-Newtonian elliptic orbit given by eq.(127) 
or (129) for e 2 < does not have a corresponding counterpart for e 2 > [see 
eq.(121)]. 

Thus Newtonian mechanics is recovered from general relativity if we ignore 
terms of order s 2 and higher for e 2 > 0, and if we ignore terms of order s and 
higher for e 2 < 0, and not from setting s = 0. If we set s = 0, then eq.(18) from 
Solution (Al) [for < e 2 < oo] gives q — oo that can be interpreted as an orbit 
of infinite radius, and eq.(35) from Solution (A2) [for < e 2 < oo] and eq.(43) 
from Solution (B) [for — oo < e 2 < 0] both give the same terminating orbit 
q = cos(0/2) for which the particle starts from the Schwarzschild horizon q = 1 
and terminates at the center of the black hole. While these orbits may not mean 
much by themselves, they are the proper limits of the corresponding orbits as 
s 2 — > 0. The small s region for the terminating orbits given by Solutions (A2) 
and (B) is non-Newtonian because the initial position of the particle is very 
close to the black hole. 

One can raise the question of whether, instead of using the energy parameter 
e and eqs.(18), (35) and (43) for the orbital trajectories, it is more convenient 
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to start with an orbit equation analogous to the Newtonian orbit equation with 
an eccentricity e appropriately defined. This is the approach used by Darwin 
and extended by Chandrasekhar that we shall discuss and compare with our 
analysis in the next section. 



9 The Analyses of Darwin and Chandrasekhar 

Chandrasekhar [5] extended Darwin's work [6] and classified the orbits from the 
solutions of eq.(8) into three types: (a) The radial geodesies (b) The bound 
orbits (k 2 < 1) and (c) The unbound orbits (k 2 > 1). The bound and un- 
bound orbits arc divided according to three categories (i), (ii) and (iii) which 
he called orbits of the first and second kind, and orbits with imaginary eccen- 
tricities, respectively. Then there are various special cases within each category. 
Chandrasekhar used various relations and inequalities for characterizing each 
type of orbit. While all possible orbits (albeit expressed in completely different 
expressions from ours) are in indeed in place, they are not clearly located on a 
"map" with coordinates that can be easily identified with the commonly used 
physical parameters. By briefly describing Chandrasekhar's work and analysis 
of these orbits and comparing them with ours, we intend to clarify this point in 
this section. 

We do not discuss the radial geodesies (of zero angular momentum) which is 
a simple special case. For the purpose of comparing Chandrasekhar's work with 
ours, it is easier to group the orbits according to the three categories (i), (ii) and 
(iii) which can be identified to those corresponding to what we call Solutions 
(Al), (A2) and (B). 

We first describe the orbits of the first kind which, as we shall see, correspond 
to the non-terminating orbits in Region I (solution Al) given by eq.(18) in 
our case. For convenience of comparison, we use our notations to describe 
Chandrasekhar's solution and parameters. To make it resemble the equation 
for the Newtonian orbit, his solution of eq.(8) is expressed as 

u= -^(1 + ecosx), (135) 

where I is the (generalized) "latus rectum" and e is the (generalized) "eccen- 
tricity" that are defined from the three real roots ui,u 2 , u 3 of the cubic equation 
under the condition A < with u\ < u 2 < u 3 as 



1 , 

u 2 = j(l +S), 

u *=4m-\ (i36) 

The "polar angle" \ m cq.(135) is related to the polar angle <j) of the spherical 
coordinate used in eqs.(3)-(5) by 
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(1 — 6/x + 2fie) 

where 



»=™> (138) 

and F(0, k) is the incomplete elliptic integral of the first kind with modulus 
k given by 

k 2 - ^\ ■ (139) 
1 - 6^ + 2^ie y ' 

The three roots of the cubic equation 

au 3 -u 2 + (2GM/h 2 )u + c 2 (k 2 - l)/h 2 = (140) 
given by eq.(136) are related to our three roots of the cubic equation (16) 

by 

au 3 = 1 - 4[i = * + 4ei, 

au 2 = 2/i(l + e) = i + 4e 2 , 
o 

qui =2/i(l-e) = ^+4e 3 , (141) 

where the dimcnsionlcss parameter /i is related to largest root e\ given in 
eq.(21) in our formulation by 

»=\- e ^ ( 142 ) 

and where the geometrical eccentricity e is related to e 2 and e 3 in our for- 
mulation by cq.(26). However, Chandrasekhar not only extended the range of 
values of e to beyond 1 but also used it for the orbits of the second kind. It 
should be noted that for orbits of the first kind, e has a simple geometrical mean- 
ing only when < e < 1. For orbits of the second kind, e has no geometrical 
meaning for any value of e. 

Indeed characterizing the orbits of the first and second kinds according to 
(fx, e) is equivalent to characterizing these orbits according to (ei, e 2 ) since e 3 = 
— (ei + e 2 ). The relationships of the parameters /i, e and our parameters e 2 , s 2 
can be shown to be given by the following equations: 

s 2 = M [l- M (3 + £ 2 )], (143) 

and 

o (1 - 4u)(l - e 2 ) 

1-/1 +£ 2 2 V ' 
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It is clear from these two equations that rather complicated combinations of 
the parameters /i and e are required to express the parameters s 2 and e 2 which 
are related simply to the angular momentum and total energy of a particle 
in a gravitational field. These dimcnsionlcss energy and field parameters e 2 
and s 2 arise naturally in the analytic solutions of eq.(13). A large part of 
Chandrasekhar's work includes expressing the constraints on [i and e that are 
required for different types of orbits. An inequality which is repeatedly stated 
that must hold for the orbits of the first and second kinds, bound or unbound, 
is 1 — ^(3 + e 2 ) > 0, which, as can be seen from eq.(143), is simply a requirement 
that s 2 > 0. Another required inequality for the orbits of the first and second 
kinds, bound or unbound, is 1 — 6/i — 2[ie > 0, which can be shown to be simply 
ei > e 2 . Other required inequalities such as fi < 1/4, < e 2 < 1 for the bound 
orbits, and /i < 1/4, e 2 > 1 for the unbound orbits, for the first and second 
kinds, are clearly associated with the identifications < e 2 < 1 and e 2 > 1 
respectively. We again mention that while e has a simple physical meaning of 
geometrical eccentricity for the orbits of the first kind in the range < e < 1, 
it does not have any simple physical meaning for the other two kinds of orbits 
for any value of e; whereas our parameter e 2 defined by eq.(15) not only has a 
simple physical meaning for all values in the range — oo < e 2 < +oo and for all 
orbits, it also has a well-defined Newtonian limit as e — > ejy given by eq.(109) 
in the limit of very small s values for the range < e < oo. Other relationships 
such as e — implies k 2 = 0, and 2/x(3 + s) = 1 implies k 2 = 1 can be seen 
easily from eq.(139). 

For the case of stable circular orbits arising from the orbits of the first kind 
for which e = 0, we have \x = (2q c )~ 1 , where 3 < q c < oo is given by eq.(58), 
and 



In terms of our parameter s' 2 , we have 



and the corresponding value of e 2 is 




(147) 



The parameter fj, can be expressed in terms of s' 2 and the corresponding e 2 




(145) 



M = 



^{l-(l-12 S f) 1/2 } 



l-3e 2 -(l + 3e 2 ) 1 / 2 
9(1 -e 2 ) 



(148) 
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For ISCO (the innermost stable circular orbit), q c = 3, [i = 1/6, k 2 = 8/9, 
e 2 = -1/3 and s 2 = 1/12. 

For the orbits of the second kind which correspond to the terminating orbits 
in Region I (Solution A2) in our case, the roots of the cubic equation (140) are 
still given by the same eq.(141), and the parameters \x and e are still defined or 
related to our formulation by eqs.(142) and (26), but the orbit equation is given 

by 

•-(5SsH) + (55If- i 7 £ ) t "4 <" 9) 

where £ is related to the polar angle (f> by 

1 ^(ik), (150) 



r (1 - 6^ + 2 y ue) 1 /2 ^2' 

where the modulus of the incomplete elliptic integral of the first kind F(0, k) 
is given by the same expression (139). Equations (143), (144), and the state- 
ments made following eq.(144) all hold for the orbits of the second kind. Re- 
member that the orbit equation (149) for the special case e = and hence 
k 2 = still gives, not a circular orbit, but a terminating orbit given by eq.(67). 
On the other hand, the orbit equation (149) for the special case 2//(3 + e) = 1 
and hence k 2 = 1 gives an unstable circular orbit with a radius given by eq.(69), 
and \x can be expressed in terms of the radius of this circle by [i = (q u — l)/ (4g„) 
or 1/4— l/(4g„), where 1.5 < q u < 3. The expressions for k 2 , s\ and e\ in terms 
of q u are remarkably similar to those in eqs.(145)-(147), and they are given by 



and 



k2 = x ^ (lgl) 

1 2q u 

si = (152) 
^qu 



el = J^- 2 , (153) 



(I - Qu) 

while the expression for fj, in terms of s\ and the corresponding e 2 is 



1 1 f / 9x1/21 1 l-3e 2 + (l + 3e 2 ) 1 / 2 , N 

^=4-l 2 { 1+ ( 1 - 12 ^ 1 = 4 18(1 -e 2 ) • ^ 

We also have the following expression for k 2 in terms of e: 

1-k 2 = ^J, (155) 
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and the expressions for s\ and the corresponding e 2 in terms of e are given 
in eqs.(91) and (89). 

As we pointed out in eqs.(64) and (65) in our discussion of Solution (Al) 
for k 2 — 1, the geometrical eccentricity e has a physical meaning expressed by 
eq.(25) for g max and g m i„ of the orbit of the first kind (Solution Al) for < e < 1 
or —1/3 < e 2 < 1. We have used the same e defined by eq.(26) for the orbit 
of Solution (A2) to express q u which is also the radius of the asymptotic circle 
given by eq.(65) or (69). In Section 6, we described the relationship between q u 
given by eq.(69) and g m i n given by eq.(65) and the continuations of the orbits 
across the boundary from Region I to II. 

For both orbits of the first and second kinds for the special case of 2^(3+e) = 
1, e 2 and s 2 are related to e defined by eq.(26) by eqs.(89) and (91) which show 
that e 2 — > oo and s 2 — > as e — > 3 and thus the values of e are restricted to 
< 3 for the special case 2/z(3 + e) = 1. The inequality given by eq.(72) implies 
< e < 3 from eq.(71). 

For the orbits with imaginary eccentricities which correspond to the termi- 
nating orbits in Region II in our case, the roots of the cubic equation (140) 
consist of one real root and two complex conjugate roots and are related to the 
roots in our formulation by 



au\ = 1 — A:[i = ^ + 4a, 
au 2 = 2/z(l + is) = ^ + 46, 

au 2 = 2/z(l - ie) = \ + 46, (156) 

where a, 6 and 6 are defined by eq.(40) and the line below it. The parameters 
H and e are related to a, 6 and 6 and to A and B of eqs.(39) and (40) by 

V= 1 --a= l --{A + B), (157) 

and 

£ - 1 l-6a " l-Q{A + By (158) 

As in the case e > 1 for Solution (Al) and for all values of e defined by 
eq.(26) for Solution (A2), the £ here defined by eq.(158) for Solution B has no 
simple physical meaning at all. The orbit equation is expressed in the form 

u= i(l + e tan|), (159) 

where £ is related to the polar angle <p by an elliptic integral, and the squared 
modulus k 2 of the elliptic integral is related to £ and \x in a somewhat more 
complicated way through several equations (see [5] p. 112). The relationships 
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of the parameters fi, e and our parameters e 2 ,s 2 are given by the following 
equations: 



.s 2 = n[l - ^(3-e 2 )], (160) 

and 



2 _ (l-4 M )(l + e 2 ) 



1 ~ e - [i-Ms-^r (161) 

The bound orbits for orbits with imaginary eccentricities are specified by 
/U < 1/4 and the unbound orbits by \x > 1/4 and 1 — /u(3 — e 2 ) > 0, and there 
is no upper limit for e 2 for either bound or unbound orbits. The parameter e is 
related to the three roots a, b, and b of the cubic equation (16) by eq.(158). The 
bound and unbound orbits are specified simply in our case by e 2 < 1 (which 
ensures that a > —1/12 for cq.(43)) and e 2 > 1 (which ensures that a < —1/12). 

The motivation of Darwin's and Chandrasckhar's analyses is to express the 
orbit equation (135) in a way which resembles the Newtonian orbit equation 
with the introduction of the "latus rectum" I (and thus fi) and "eccentricity" 
e. However, this is done at the expense that not only the "polar angle" \ or 
£ is not easily related to the polar angle <j) 01 the spherical coordinate, but 
also that I and e (that can become imaginary) do not have simple physical 
meanings for most of these orbits. The two parameters t and e in essence 
represent roots of the cubic equation (140), and the entirety of all possible orbits 
cannot be conveniently placed in a parameter space that is directly related to 
the commonly used physical quantities. 

Incidentally, for the unbound orbits of the first kind (Solution Al for e 2 > 1) 
and the unbound terminating orbits in Region II (solution B for e 2 > 1), an 
impact parameter b can be defined. Chandrasekhar gave this impact parameter 
incorrectly as (see [5], p. 114) 

b= i^f_ . (162) 
V« 2 - 1 V ' 

The correct impact parameter is [3] 

&=-#T> ( 163 ) 
Vr - 1 

where h and k are defined by eqs.(5) and (6). In terms of e 2 and s 2 , the 
dimensionless impact parameter is given by eq.(34). The possible cause of error 
in eq.(162) is the following oversight: When a particle of rest mass m initially 
very far away from the star or black hole has an initial velocity and an 
impact parameter b from the star or black hole, the angular momentum h per 
unit rest mass of the particle is not h — v^b, but should be h = mvoob/mo, 
where m = m a / yjl — v^/c 2 . 
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10 Summary and Commentary 



We have presented three possible analytic solutions for the orbit of a particle 
in the Schwarzschild geometry in explicitly simple forms by eqs.(18), (35) and 
(43), and we have presented a universal map with coordinates (e 2 ,s) [as in [l] 
and [2] on which all possible orbits arising from these three solutions can be 
placed, where the two dimensionless parameters e 2 and s 2 [defined by eqs.(15) 
and (11)] which we call the energy and field parameters respectively, are in 
the ranges —00 < e 2 < 00 and < s < 00. We have presented but not used 
two other maps with coordinates (n 2 ,s 2 ) (13) and (—53,-52) (14) that could 
be interesting alternatives where k, 172 and 33 are defined by eqs.(6) and (14). 
The non-terminating orbit solution (18) and the terminating orbit solution (35) 
occur in Region I (A < 0), and the other terminating orbit solution (43) occurs 
in Region II (A > 0), where A is defined by eq.(17). In contrast to the analyses 
of other authors that used various mathematical inequalities for separating and 
discussing the different regions, Region I is explicitly shown (in[T]) to be bounded 
on the left by s — s[ given by eq.(48) that is obtained from setting A = and 
k 2 = [see eqs.(20) and (42)], and it is bounded on the top by s = si given 
by eq.(50) that is obtained from setting A = and k 2 = 1, and Region II is 
bounded on the top for e 2 < 1 by s = S2 given by eq.(52). 

The use of e 2 and s 2 for characterizing an orbit is seen to arise naturally from 
eq.(10) or (13), and the advantage of this characterization is enhanced by the 
fact that the non-terminating orbits inside Region I (excluding its boundaries) 
can be conveniently classified into three types: the elliptic-type for —1/3 < 
e 2 < 1 J3]), the parabolic-type for e 2 = 1 Q, and the hyperbolic-type for 
e 2 > 1 With the star or black hole at the origin, the elliptic-type orbit 
is a bound precessing elliptic orbit with the precessional angle A</> given by 
eq.(27) and with (j max and g m i n given by eqs.(23) and (24). The parabolic-type 
orbit is an unbound orbit in which a particle approaches from infinity at (f> = 
and can, in general, go around the star or black hole a number of times as 
its polar angle increases from to 2K(k)/'y before going off to infinity. The 
hyperbolic-type orbit is also an unbound orbit in which a particle approaches 
from infinity along an incoming asymptote at an angle "Ji given by eq.(29) to 
the horizontal axis, turning counter-clockwise about the star to its right on the 
horizontal axis, and leaving along an outgoing asymptote at an angle \&2 given by 
eq.(30). The x coordinates of the points where the two asymptotes intersect the 
horizontal axis relative to the origin are given by eq.(32), a straight line through 
the origin that makes an angle x with the horizontal axis given by eq.(33) is 
the symmetry axis of the hyperbolic-type orbit, and the impact parameter or 
the perpendicular distance from the origin to the asymptote is given by eq.(34). 
As for the terminating orbits inside Region I, a particle initially at a distance 
qi given by eq.(36) from the black hole at <fi — plunges into the center of the 
black hole when the polar angle <j) = <j)\ given by eq.(37), and q\ is finite for all 
values of e 2 (JgJ [7] and |8| . Region II allows only terminating orbits. The initial 
distance q2 of the particle from the black hole is finite and is given by eq.(44) 
for — 00 < e 2 < 1 ([9] and 12), it is infinite at (f> = for e 2 = 1 (10), and it is 
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infinite at tfi — ^2 given by eq.(46) for e 2 > 1 (11). The particle plunges into 
the center of the black hole when its polar angle (f> = 4>2 given by eq.(45). 

Two types of orbits occur on the boundary : stable circular orbits of radii 
q c given by eq.(57) and terminating orbits given by eq.(67). Two types of orbits 
occur on the boundary Sj: unstable circular orbits of radii q u given by eq.(69), 
and asymptotic orbits given by eq.(62) with the particle starting from q — ^ max 
for —1/3 < e 2 < 1 (15 1, or with the particle starting from infinity at the polar 



angle <j) = \&i given by eq.(66) for e > 1 (16 and 17). The asymptotic orbits end 



up circling the star or black hole with a radius that asymptotically approaches 
9min given by eq.(65). As the boundary s\ is crossed from Region I to Region 
II, the non-terminating orbit of Region I becomes an asymptotic orbit on s% 
and then becomes the terminating orbit of Region II, while the terminating 
orbit of Region I becomes an unstable circular orbit on s\ and becomes part of 
the terminating orbit in Region II. On the other hand, as the boundary s[ is 
crossed from Region II to Region I, the terminating orbit of Region II becomes 
a terminating orbit on and then becomes the terminating orbit of Region 
I, while a stable circular orbit becomes possible on and then becomes an 
elliptical orbit in Region I. These transitions are illustrated by pictures of the 
corresponding changes in the "effective potential" curve W ([I8j23|. All these 
observations could not have been so clearly presented without the help of the 
two curves that clearly mark the boundaries s' x and Si between Regions I and II 
on a universal map ([I] and [2]) . We consider the presentation of all these clearly 
identifiable transitions of orbits from Region I to Region II to be one of the 
distinctive results of this paper. 

The use of e 2 for characterizing the orbits is also related to the fact that in 
the weak field limit and for e 2 > 0, e becomes the Newtonian eccentricity e^r 
that describes the classical Newtonian orbit given by eq.(llO) or (121), as can be 
seen from eqs.(107)-(109) and (112). Thus the usual division of the Newtonian 
orbit into elliptic for < ejv < 1, parabolic for ejy = 1, and hyperbolic for 

> 1 is seen to be a special case of the division for the general relativistic 
non-terminating orbits in Region I. We showed that in Region I close to s[, 
there exist elliptic orbits with a non-Newtonian eccentricity as shown in eqs. 
(127) and (129). For the hyperbolic-type orbits in the Newtonian limit, an 
approximate analytic expression [eq.(134)] has been presented that gives the 
difference between the general relativistic and the Newtonian calculations for 
the angle between the two asymptotes and this expression is different from the 
one given by other authors [13]. 

The power of our analytic expressions for the orbits in the Schwarzschild 
geometry is further illustrated by the numerous analytic simple expressions for 
many quantities that are presented for various special cases in Section 7. These 
examples include the cases of k 2 = 0, 1/2, and 1, and of e 2 = 1 and —1/3, 
and the expressions include those for e, the geometrical eccentricity defined 
by eq.(25), qi and 92 discussed above, and also 92 max in the region e 2 < 0. 
Every expression from Eq.(73) to Eq.(106) has been used to check the extensive 
numerical data for various orbits that were presented in refs.[2] and [3], and we 
consider these simple relationships to be a bonus of our analytic formulation of 
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the problem. 

To better appreciate our presentation and characterization of all possible 
orbits in the Schwarzschild geometry, we described in Section 9 the analysis of 
Chandrasekhar [5] in terms of our notations for the purpose of comparisons. 
Equations (135), (149) and (159) are, respectively, mathematically equivalent 
to our eqs.(18), (35) and (43). Instead of using the parameters e 2 and s as we 
adopted, Chandrasekhar's formulation essentially used the parameters e defined 
by eqs.(26) and (158), and // defined by eq.(138), for characterizing the orbits, 
but he often had to refer to the total energy [E of eq.(6)] and total angular 
momentum L [which is a constant multiple of h of eq.(5)] of the particle and 
the effective potential energy curve when discussing certain features of the or- 
bits. Although the orbit equations appear to involve only simple trigonometric 
functions, the "polar angles" x an d £ are not the actual physical polar angle </> 
but are complicated implicit functions of cj> [see eqs.(137), (150), and the quote 
after eq.(159)]. The two important physical quantities, E and h, are not only 
complicated functions of e and /i, but the relations are different for different 
types of orbits [see eqs.(143), (144), (160), and (161)]. The parameter e can be 
real or purely imaginary, but it has a physical meaning given by eq.(25) only 
in the range of real values between and 1. The boundaries between different 
regions for different types of orbits are expressed as certain equations that relate 
e and /i. Chandrasekhar's formalism does not present all possible orbits on a 
single universal map. The numerous equations in Section 9, however, provide 
the connections of Chandrasekhar's formulation to ours that would be helpful 
to many readers who are familiar with the former. Equation (163) or (34) also 
corrects an erroneous expression (162) given by Chandrasekhar. 
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Fig. 1 The boundaries of Regions I, II, and II' are displayed. The vertex 
point V occurs at the intersection of the left boun dary and the upper 
boundary s^. Circular orbits are represented on these two boundaries of 
Region I with the innermost stable circular orbit occurring at V. Unbound 
orbits are possible when e 2 is greater than or equal to 1 . 



Figure 1: The boundaries of Regions I, II and IF are displayed. The vertex point 
V occurs at the intersection of the left boundary s' x and the upper boundary si. 
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Fig. 2 Curves of constant Within Region I only three values ofli? are 
displayed. All curves intersect the point V. The dashed horizontal line 
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Figure 2: Curves of constant k 2 . Within Region I only three values of k 2 are 
displayed. All curves intersect the point V. The dashed horizontal line through 
V has s 2 = 1/12. 
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Fig. 3 Elliptic-type orbit in Region I for = 0.25 and s = 0.2 



gure 3: Elliptic-type orbit in Region I for e 2 = 0.25 and s = 0.2. 
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Fig. 4 Parabolic-type orbit in Region I for e~ = 1 and 5 = 0.2 . 



Figure 4: Parabolic-type orbit in Region I for e 2 = 1 and s = 0.2. 
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Fig. 5 Hyperbolic orbit in Region I fore 2 =4 and,s = 0.2 . The 
asymptotes to the orbit are also displayed. 



Figure 5: Hyperbolic-type orbit in Region I for e 2 = 4 and s = 0.2. The 
asymptotes to the orbit are also displayed. 
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Fig. 6 Terminating orbit in Region I for e 2 = 0.25 and s = 0.2. 



;ure 6: Terminating orbit in Region I for e 2 = 0.25 and s = 0.2. 

50 




51 




52 




Fig. 9 Terminating orbit in Region II for ^ = 0.25 and s = 0.3 . 



Figure 9: Terminating orbit in Region II for e 2 = 0.25 and s = 0.3. 
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Figure 11: Terminating orbit in Region II for e 2 = 4 and s = 0.3. The asymptote 
to the orbit is also displayed. 



■55 




5G 




Fig. 13 (k 2 , s 2 ) parameter space. The thick solid lines marked by A 2 = 
and A 2 = 1 represent the boundaries of Region I and the dashed solid lines 
from left to right represent curves of constant e 2 with values of -10, -5, 
-2, 0, 1, 2, 5, and 10 respectively with the highlighted vertical dashed line 
representing e 2 = 1 . 



Figure 13: (k, 2 ,s 2 ) parameter space. The thick lines marked by k 2 = and 
k 2 = 1 represent the boundaries of Region I and the dashed solid lines represent 
curves of constant e 2 with values of, from left to right, —10, —5, —2, 0, 1, 2, 5, and 
10 respectively with the highlighted vertical dashed line representing e — 1. 
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Fig. 14 Curves of constant e 2 using g 2 andg- g as coordinates. The thick 
symmetric curves below the horizontal axis are the boundaries of Region 
I and represent A 2 = 1 to the right of the vertical axis and k 1 = to the left 
of the vertical axis. The thick dashed line represents the boundary of 
Region II. 



Figure 14: Curves of constant e 2 using gi and 53 as coordinates. The thick 
symmetric curves below the horizontal axis are the boundaries of Region I and 
represent k 2 — 1 to the right of the vertical axis and fc 2 = to the left of the 
vertical axis. The thick dashed line represents the boundary of Region II. 
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Figure 16: Asymptotic orbit for e 2 = 1 and k 2 — 1 and s\ = 0.25. 
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Fig. 17 Asymptotic orbit for e 2 = 4 and A 2 = 1 and s 1 = 0.2210354 . 
The asymptote to the orbit is also displayed. 



Figure 17: Asymptotic orbit for e 2 = 4 and k 2 = 1 and si = 0.2210354. The 
asymptote to the orbit is also displayed. 
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;ure 19: W curve for e 2 = 0.25 and s = 0.2648124. 
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Figure 20: W curve for e 2 = 0.25 and s = 0.3. 
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